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FOREWARD 


The  NRL  Plasma  Formulary  originated  over  twenty  five  years  ago  and 
has  been  revised  several  times  during  this  period.  The  guiding  spirit  and  per¬ 
son  primarily  responsible  for  its  existence  is  Dr.  David  Book.  I  am  indebted  to 
Dave  for  providing  me  with  the  TgX  files  for  the  Formulary  and  his  continued 
suggestions  for  improvement.  The  Formulary  has  been  set  in  TgX  by  Dave 
Book,  Todd  Brun,  and  Robert  Scott.  Finally,  I  thank  readers  for  communicat¬ 
ing  typographical  errors  to  me. 


2 


CONTENTS 


Numerical  and  Algebraic  . 4 

Vector  Identities . 5 

Differential  Operators  in  Curvilinear  Coordinates . 7 

Dimensions  and  Units . 11 

International  System  (SI)  Nomenclature . 14 

Metric  Prefixes . 14 

Physical  Constants  (SI) . 15 

Physical  Constants  (cgs)  . 17 

Formula  Conversion  . 19 

Maxwell’s  Equations  . 20 

Electricity  and  Magnetism . 21 

Electromagnetic  Frequency/ Wavelength  Bands . 22 

AC  Circuits . 23 

Dimensionless  Numbers  of  Fluid  Mechanics  . 24 

Shocks  . 27 

Fundamental  Plasma  Parameters . 29 

Plasma  Dispersion  Function . 31 

Collisions  and  Transport  . 32 

Ionospheric  Parameters . 41 

Solar  Physics  Parameters  . 42 

Thermonuclear  Fusion  . 43 

Relativistic  Electron  Beams . 45 

Beam  Instabilities  . 47 

Approximate  Magnitudes  in  Some  Typical  Plasmas . 49 

Lasers  . 51 

Atomic  Physics  and  Radiation . 53 

Atomic  Spectroscopy . 59 

Complex  (Dusty)  Plasmas . 62 

References . 66 


3 


NUMERICAL  AND  ALGEBRAIC 

Gain  in  decibels  of  P2  relative  to  P± 

G  =  101og10(P2/Pi). 


To  within  two  percent 

(2tt)1/2  «  2.5;  7i2  «  10;  e3  «  20;  210  «  103. 


Euler-Mascheroni  constant1  7  =  0.57722 


Gamma  Function  T(x  +  1)  =  xF(x): 


r(l/6)  =  5.5663  r(3/5) 

r(l/5)  =  4.5908  T(2/3) 

r(l/4)  =  3.6256  r(3/4) 

r(l/3)  =  2.6789  r(4/5) 

r(2/5)  =  2.2182  r(5/6) 

r(l/2)  =  1.7725  =  r(l) 


1.4892 

1.3541 

1.2254 

1.1642 

1.1288 

1.0 


Binomial  Theorem  (good  for  |  x  |<  1  or  a  =  positive  integer): 


( l+x)a  = 


E  Q)x" =i+a*+ 


a(a  —  1)  2  ol{ol  —  l)(a 

■x  -\ - 


2! 


3! 


2)  3  , 

—  x  +  .  .  . 


k  =  0 


Rothe-Hagen  identity2  (good  for  all  complex  x,  y,  z  except  when  singular): 


Ex  /x  +  kz\ 
x  +  kz  V  k  ' 


y 


x  -(-  kz  '  k  '  y  +  (n  —  k)z 


k= 0 


/■y  +  (n  —  k)z' 

^  n  —  k  2 

x  +  y  (X  +  y  +  nz\ 


x  +  y  +  nz 


+  +  nz\ 

'  n  ' 


Newberger’s  summation  formula3  [good  for  /i  nonintegral,  Re  (a  +  /?)  >  —  1] 

OO 

(  1)  Jot  —  -yn  (■2;)  J(3-\-/yn  (■2') 


E 


n  =  —  00 


7T 


n  +  /i 


Sin  [ITT 


'  Jat+'y/j,  (z)Jf3  —  ^fi  (z)- 
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VECTOR  IDENTITIES4 

Notation:  /,  g ,  are  scalars;  A,  B,  etc.,  are  vectors;  T  is  a  tensor;  /  is  the  unit 
dyad. 

(1) A-BxC  =  AxB-C  =  B-CxA  =  BxC-A  =  C-AxB  =  CxA-B 

(2)  A  x  (B  x  C)  =  (C  x  B)  x  A  =  (A  •  C)B  -  (A  •  B)C 

(3)  A  x  (B  x  C)  +  B  x  (C  x  A)  +  C  x  (A  x  B)  =  0 

(4)  (A  x  B)  ■  (C  x  D)  =  (A  •  C)(B  •  D)  -  (A  •  D)(B  •  C) 

(5)  (A  x  B)  x  (C  x  D)  =  (A  x  B  •  D)C  -  (A  x  B  •  C)D 

(6)  V (f g)  =  V (gf)  =  fVg  +  gVf 

(7)  V  •  (/A)  =  /V  •  A  +  A  •  V/ 

(8)  V  x  (/A)  =  /V  x  A  +  V/  x  A 

(9)  V  •  (A  x  B)  =  B  •  V  x  A  —  A  •  V  xB 

(10)  V  x  (A  x  B)  =  A(V  •  B)  -  B( V  •  A)  +  (B  •  V)A  -  (A  •  V)B 

(11)  A  x  (V  x  B)  =  (VB)  ■  A  -  (A  •  V)B 

(12)  V (A  •  B)  =  A  x  (V  x  B)  +  B  x  (V  x  A)  +  (A  •  V)B  +  (B  •  V) A 

(13)  V2/  =  V  •  V/ 

(14)  V2A  =  V(V  •  A)  -  V  x  V  x  A 

(15)  V  x  V/  =  0 

(16)  V  •  V  x  A  =  0 

If  ei,  e2,  e3  are  orthonormal  unit  vectors,  a  second-order  tensor  T  can  be 
written  in  the  dyadic  form 

(17)  T  =  ^  Tij  Qi  Gj 

In  cartesian  coordinates  the  divergence  of  a  tensor  is  a  vector  with  components 

(18)  (V-  T)i  =  E.tSTWari) 

[This  definition  is  required  for  consistency  with  Eq.  (29)].  In  general 

(19)  V  •  (AB)  =  (V  •  A)B  +  (A  •  V)B 

(20)  V  •  (/T)  =  V/-T+/V-T 
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Let  r  =  ix  +  jy  +  k z  be  the  radius  vector  of  magnitude  r,  from  the  origin  to 
the  point  x,y,z.  Then 


(21) 

V  •  r  =  3 

(22) 

V  X  r  =  0 

(23) 

Vr  =  r/r 

(24) 

V(l/r)  = 

—  r/r3 

(25) 

<1 

i-T 

CO 

=  47r5(r) 

(26) 

Vr  =  / 

If  V  is  a  volume  enclosed  by  a  surface  S  and  dS  =  n dS ,  where  n  is  the  unit 
normal  outward  from  V, 


(27)  /  dVVf=  /  dSf 

J  v  J  s 

(28)  /  dV  V  •  A  =  /  dS  ■  A 

J  v  J  s 

(29)  /  dVV  T  =  /  dS  ■  T 

J  v  J  s 

(30)  /  dV V  x  A  =  /  dS  x  A 

J  v  J  s 

(31)  /  dV{f\/2g  -  g\72 f)  =  /  dS-(fVg-gVf) 

J  v  J  s 

(32)  [  dV{  A  -  VxVxB  —  B-VxVx  A) 


v 


=  /  dS  ■  (B  x  V  x  A  -  AxVxB) 


If  S1  is  an  open  surface  bounded  by  the  contour  C ,  of  which  the  line  element  is 

dl, 


(33)  j  dS  x  V/  =  (t)  dlf 
s  J  c 
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(34) 

(35) 

(36) 


dS  •  V  X  A  =  (7)  dl  ■  A 

s  J  c 

(dS  xV)xA  =  (j)  dl  x  A 
s  J  c 


dS  ■  (V/  xVg)=  (h  fdg 
s  J  c 


gdf 


DIFFERENTIAL  OPERATORS  IN 
CURVILINEAR  COORDINATES5 


Cylindrical  Coordinates 


Divergence 


V  •  A 


Id  1  dA^  dAz 

- (rAr)  - £  H - 

r  dr  r  dcf)  dz 


Gradient 


(V/)r 


2i 

dr  ’ 


(V/)* 


1 21 

r  dcj)  ’ 


(V/), 


21 

dz 


Curl 


(V  x  A)r 


1  dAz 

r  dcj) 


dA ^ 
dz 


(V  x  A)0 


dAr  dAz 
dz  dr 


(V  x  A)z 


-|-(rA*) 


1  dAr 
r  dcj) 


Laplacian 


V2/ 


I A  ,  dlL 

r  dr  \  dr  y  r2  dcj)2  dz2 
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Laplacian  of  a  vector 


o  2  2  dA(k 

(V2A  )r  =  v2^-— — f 

r21  o</>  r21 


2  s  2  2  5Ar 

(v2a)^,  =  V2A^  +  —  — f. 

r22  dcp 


(V2A),  =  v2a. 


Components  of  (A  ■  V)B 


5Sr 

A<^ 

5Br 

dr 

1 

r 

dcf) 

OB $ 

_1_ 

A^ 

dB $ 

dr 

1 

r 

dcf) 

dBz 

_l_ 

A^ 

dBz 

dr 

1 

r 

dcf) 

5.Br 


.  dB^  Br 

'  Z  r\  ' 

oz  r 


Divergence  of  a  tensor 

(V  ■  T)r  =  -^-(rTrr)  +~^  +  ^  ~  ^ 

r  or  r  ocp  oz  r 


rx-7  T\  _  ^  (  rr  \  \  |  dTzcf,  |  T<£r 

(V  ■  /  )<£  —  -  —\rTr^)  -\ - — - 1 - - - 1 - 

r  Or  r  Ocp  Oz  r 


15  1  dTsz  dTzz 

(V  .  T)*  =  -  — (rT„)+  ^  +  — - 

r  or  r  ocp  oz 
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Spherical  Coordinates 

Divergence 


V  •  A 


J__2_(r2Ar)  + 

r2  dr  r  sin  0  80 


1  d  , 

—  (sm  6Ae)  + 


1  dA# 

rsind  dcj) 


Gradient 


(V/)r 


(V/)„ 


19/. 

r  d6  ’ 


(V/)* 


1  9/ 


r  sin  d  d0 


Curl 


(V  x  A)r 


i  a 

——  (sm  OAfi) 


r  sin  0  dO 


1  dAe 

r  sin  0  dcf) 


(V  x  A)e 


1  dAr 
rsind  dcf) 


(V  x  A)^ 


Id  1  <9Ar 

-  — (rA0) - — 

r  dr  r  dO 


Laplacian 


J_L  (r2A)  +  1  8  +  _ 1 _ 

r2  dr  \  dr  J  r2  sin  0  dO  \  dO  J  r2  sin2  0  dcj)2 


Laplacian  of  a  vector 


(V2  A)r 


V2  Ar 


2  dAe 
r2  dO 


2  cot  0  A 
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dA 


4> 


r2  sin  0  dcj) 


2  2  2  dAr 

(V2A)ff  =  V2A0  +  — 


Af 


2  cos  0  dA 


<t> 


r2  dO  r2  sin2  0  r2  sin2  0  dcj) 


(V2A)0  =  v2a^ 


Arf,  2  dAr  2  cos  0  dAe 

~ t-  o  n  i 


r2  sin2  0  r2  sin  0  dcj)  r2  sin2  0  dcj> 
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Components  of  (A  •  V)B 


,  a  wt=>\  a  &Br  Aq  8Br  dBr  AqBq  +  A^B^ 

(A  *  VJjJr*  —  J±r  “  i  “ ”  i 


dr  r  89  r  sin  9  dcf) 


r 


/a  WON  a  dBe  ,  Ae  dBe  A ^  8B0  AeBr  cot  9A<t>B<j> 

(A-VB  )&  =  Ar— - 1- - —  +  — r^-^  + 


dr  r  89  rsinf?  dcf> 


r 


r 


/a  a  dB(t>  A&  dB  ^  A ^  8B ^  A^Br  cot  9Art>B0 

(A  •  VB)<£  =  Ar— - 1 - — 1 - ; — -  — — 1 - r 


dr  r  89  rsin0  dcf> 


r 


r 


Divergence  of  a  tensor 

(V  ■  T)r  =  (r2Trr)  H - -  -^(sin0T0r) 

rz  dr  r  sin  9  89 


1  dT^r  Tqo  +  T<t>ct> 


r  sin  9  dcf) 


r 


(V  .  T)e 


Id.  2  v 
—  —{r2Tre)  + 

jLj  If* 


8 


r  sin  9  89 


—  (sin  9Too ) 


1  dT^g  Tor  cot  9T(p(p 
r  sin  9  dcf)  r  r 


(V  .  T)4, 


1  d  2  \  1  d 

_ (rTrh)  H _ 

r2  dr  1  r  sin  9  89 


—  (sin  9Td(f, ) 


1 


r  sin  9  dcf) 


dT T^r  cot  9T^q 


r 


r 
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DIMENSIONS  AND  UNITS 


To  get  the  value  of  a  quantity  in  Gaussian  units,  multiply  the  value  ex¬ 
pressed  in  SI  units  by  the  conversion  factor.  Multiples  of  3  in  the  conversion 
factors  result  from  approximating  the  speed  of  light  c  =  2.9979  X  1010  cm/sec 
^  3  X  1010  cm/sec. 


Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

Capacitance 

C 

tV 

ml2 

Q 

l 

farad 

9  X  1011 

cm 

Charge 

Q 

ml/2^3/2 

coulomb 

3  X  109 

statcoulomb 

t 

Charge 

density 

Q 

m1'2 

coulomb 

/m3 

3  X  103 

statcoulomb 

/cm3 

P 

l 3 

P/2t 

Conductance 

tq2 
ml 2 

l 

t 

siemens 

9  X  1011 

cm/sec 

Conductivity 

a 

tq2 
ml 3 

1 

t 

siemens 

/m 

ampere 

9  X  109 

sec  - 1 

Current 

hi 

Q 

m 1'2l3'2 

3  X  109 

statampere 

t 

Q 

t2 

m1/2 

Current 

JJ 

ampere 

3  x  105 

statampere 

l2t 

n/H2 

density 

/m2 

/cm2 

Density 

P 

m 

¥ 

m 

¥ 

kg/ m3 

10-3 

g/ cm3 

<7 

m1/2 

Displacement 

D 

coulomb 

/m2 

127 T  x  105 

statcoulomb 

/cm2 

l2 

l^2t 

Electric  field 

E 

ml 

m1/2 

volt  /  m 

CO  1 

X 

o 

1 

st  at  volt /cm 

t2q 

l^2t 

Electro- 

motance 

Emf 

ml2 

volt 

1  _  2 

-  x  10 

3 

st  at  vo  It 

t2q 

t 

Energy 

U,  W 

ml2 

t 2 

ml2 

t 2 

joule 

107 

erg 

Energy 

density 

w,  e 

m 

T¥ 

m 

T¥ 

joule/m3 

10 

erg/ cm3 
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Dimensions 

Physical 

Sym- 

SI 

Conversion 

Gaussian 

Quantity 

bol 

SI 

Gaussian 

Units 

Factor 

Units 

ml 

ml 

Force 

F 

l2 

l2 

newton 

105 

dyne 

1 

1 

Frequency 

f,  " 

t 

t 

hertz 

1 

hertz 

ml 2 

t 

1  li 

sec/cm 

Impedance 

Z 

tq 2 

1 

ohm 

-  X  10  11 

9 

ml2 

t2 

1  li 

Inductance 

L 

q 2 

l 

henry 

-  X  10  11 

9 

sec  /cm 

Length 

l 

l 

l 

meter  (m) 

102 

centimeter 

(cm) 

1/2 

Magnetic 

H 

Q 

m  1 

ampere— 

4t t  X  10-3 

oersted 

It 

I'/H 

intensity 

turn/m 

Magnetic  flux 

4> 

ml2 

weber 

h-1 

o 

00 

maxwell 

tq 

t 

m1/2 

Magnetic 

B 

tesla 

104 

gauss 

tq 

n/2t 

induction 

Magnetic 

m,  n 

l2q 

m1/2/5/2 

O 

CO 

oersted— 

t 

t 

cXIll  L/LI  L  XiX 

moment 

cm3 

1/2 

Magnetization 

M 

Q 

m  ' 

47 T  x  10-3 

oersted 

It 

n/2t 

cXUI  L/LI  L 

turn/m 

Magneto- 

M, 

Q 

47 r 

gilbert 

t 

t2 

cXUI  L/LI  L 

10 

motance 

Mmf 

turn 

Mass 

m,  M 

m 

m 

kilogram 

O 

CO 

gram  (g) 

(kg) 

ml 

ml 

kg— m  /  s 

g— cm/sec 

Momentum 

P,P 

t 

t 

105 

Momentum 

m 

Jh 

m 

Jh 

kg/m2— s 

o 

1 

g/cm2— sec 

density 

ml 

henry/m 

1  7 

Permeability 

h 

q2 

i 

—  x  10' 

4tt 

12 


Physical 

Quantity 


Dimensions 


SI 

Units 


Conversion 

Factor 


Gaussian 

Units 


Permittivity 

Polarization 

Potential 

Power 

Power 

density 

Pressure 

Reluctance 

Resistance 

Resistivity 

Thermal  con¬ 
ductivity 

Time 

Vector 

potential 

Velocity 

Viscosity 

Vorticity 

Work 


SI 

Gaussian 

tV 

1 

ml 3 

Q 

m1/2 

l 2 

U/2t 

ml 2 

m1,2l2/2 

t2q 

t 

ml2 

ml2 

i3 

t3 

m 

m 

It 3 

It 3 

m 

m 

It 2 

It 2 

q2 

1 

ml 2 

l 

ml2 

t 

tq2 

l 

ml 3 

t 

tq 2 

ml 

ml 

T* 

"t3" 

t 

t 

ml 

mV2;1/2 

tq 

t 

l 

l 

t 

t 

m 

m 

It 

It 

1 

1 

t 

t 

ml2 

ml2 

t 2 

t2 

farad/m 

36t t  X  109 

— 

coulomb/m2 

3  X  105 

statcoulomb 

/cm2 

volt 

1  -2 
-  x  10  2 

3 

st  at  vo  It 

watt 

107 

erg/ sec 

watt/m3 

10 

erg/ cm3-sec 

pascal 

10 

dyne/cm2 

ampere— turn 
/weber 

47r  x  10~:i 

cm-1 

ohm 

-  x  10-11 

9 

sec/cm 

ohm— m 

-  x  10-9 

9 

sec 

watt  /  m— 
deg  (K) 

105 

erg/ cm— sec— 
deg  (K) 

second  (s) 

1 

second  (sec) 

weber/m 

106 

gauss— cm 

m/s 

o 

to 

cm/sec 

kg  /  m— s 

10 

poise 

s-1 

1 

sec-1 

joule 

107 

erg 
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INTERNATIONAL  SYSTEM  (SI)  NOMENCLATURE6 


Physical 

Name 

Symbol 

Physical 

Name 

Symbol 

Quantity 

of  Unit 

for  Unit 

Quantity 

of  Unit 

for  Unit 

*length 

meter 

m 

electric 

potential 

volt 

V 

*mass 

kilogram 

kg 

electric 

ohm 

Q 

*time 

second 

s 

resistance 

*current 

ampere 

A 

electric 

conductance 

siemens 

S 

*temperature 

kelvin 

K 

electric 

farad 

F 

^amount  of 

mole 

mol 

capacitance 

substance 

magnetic  flux 

weber 

Wb 

*luminous 

candela 

cd 

intensity 

magnetic 

inductance 

henry 

H 

fplane  angle 

radian 

rad 

magnetic 

tesla 

T 

jsolid  angle 

steradian 

sr 

intensity 

frequency 

hertz 

Hz 

luminous  flux 

lumen 

lm 

energy 

joule 

J 

illuminance 

lux 

lx 

force 

newton 

N 

activity  (of  a 
radioactive 

becquerel 

Bq 

pressure 

pascal 

Pa 

source) 

power 

watt 

W 

absorbed  dose 
(of  ionizing 

gray 

Gy 

electric  charge 

coulomb 

C 

radiation) 

*SI  base  unit  f Supplementary  unit 


METRIC  PREFIXES 


Multiple 

Prefix 

Symbol 

Multiple 

Prefix 

Symbol 

o 

1 

M 

deci 

d 

10 

deca 

da 

o 

1 

to 

centi 

c 

102 

hecto 

h 

o 

1 

CO 

milli 

m 

o 

CO 

kilo 

k 

10-6 

micro 

h 

106 

mega 

M 

O 

1 

CD 

nano 

n 

109 

giga 

G 

10-12 

pico 

P 

1012 

tera 

T 

10-15 

femto 

f 

1015 

peta 

P 

IQ-18 

atto 

a 

1018 

exa 

E 

PHYSICAL  CONSTANTS  (SI)7 


Physical  Quantity 

Boltzmann  constant 
Elementary  charge 
Electron  mass 
Proton  mass 
Gravitational  constant 
Planck  constant 

Speed  of  light  in  vacuum 

Permittivity  of 
free  space 

Permeability  of 
free  space 

Proton/electron  mass 
ratio 

Electron  charge/mass 
ratio 

Rydberg  constant 

Bohr  radius 
Atomic  cross  section 
Classical  electron  radius 
Thomson  cross  section 

Compton  wavelength  of 
electron 

Fine-structure  constant 

First  radiation  constant 

Second  radiation 
constant 

Stefan-Boltzmann 

constant 


Symbol 

Value 

Units 

k 

1.3807  X  10-23 

JK-1 

e 

1.6022  x  10-19 

C 

me 

9.1094  x  10-31 

kg 

mp 

1.6726  X  10-27 

kg 

G 

6.6726  X  10-11 

m3s_2kg_  1 

h 

6.6261  X  10-34 

J  s 

55+ 

5- 

to 

^1 

1.0546  X  10“34 

J  s 

c 

2.9979  X  108 

m  s“  1 

eo 

8.8542  x  10-12 

F  m_1 

Ho 

4tt  x  10-7 

Hm-1 

mp/me 

1.8362  x  103 

e/me 

1.7588  x  1011 

Ckg-1 

4 

me 

1.0974  x  10' 

m 

8eo 2ch3 

ao  =  eoh2  /  nme2 

5.2918  x  10-11 

m 

nao2 

8.7974  x  10-21 

m2 

re  =  e  /47reo me 

2.8179  x  10-15 

m 

(8tt/3  )re2 

6.6525  X  10-29 

m2 

h/mec 

2.4263  X  10-12 

m 

h/mec 

3.8616  X  10-13 

m 

a  =  e2 /2eohc 

7.2974  x  10-3 

a-1 

137.04 

ci  =  2irhc2 

3.7418  x  10“16 

Wm2 

C2  =  he/  k 

1.4388  X  10“2 

m  K 

cr 

5.6705  X  10-8 

W  m-2K-4 

Physical  Quantity 


Symbol 


Value 


Units 


Wavelength  associated 

Ao  =  he/ e 

1.2398  X  10“6 

m 

with  1  eV 

Frequency  associated 

vo  =  e/h 

2.4180  X  1014 

Hz 

with  1  eV 

Wave  number  associated 

ko  =  e/hc 

8.0655  X  105 

m_  1 

with  1  eV 

Energy  associated  with 

hv  o 

1.6022  x  10-19 

J 

1  eV 

Energy  associated  with 

he 

1.9864  x  10-25 

J 

1  m  -  1 

Energy  associated  with 

me3 /8eo2  h2 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

k/e 

8.6174  x  10“5 

eV 

1  Kelvin 

Temperature  associated 

e/k 

1.1604  x  104 

K 

with  1  eV 

Avogadro  number 

Na 

6.0221  x  1023 

mol-1 

Faraday  constant 

F  =  NAe 

9.6485  X  104 

C  mol-1 

Gas  constant 

R  =  N  Ak 

8.3145 

J  K  — 1mol  — 1 

Loschmidt’s  number 

n0 

2.6868  X  1025 

m_J 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6605  X  10-27 

kg 

Standard  temperature 

To 

273.15 

K 

Atmospheric  pressure 

po  =  n0kTo 

1.0133  X  105 

Pa 

Pressure  of  1  mm  Hg 

1.3332  x  102 

Pa 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTo/po 

2.2414  x  10“2 

m3 

Molar  weight  of  air 

Mair 

2.8971  X  10“2 

kg 

calorie  (cal) 

4.1868 

J 

Gravitational 

9 

9.8067 

ms-2 

acceleration 
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PHYSICAL  CONSTANTS  (cgs)7 


Physical  Quantity 

Symbol 

Value 

Units 

Boltzmann  constant 

k 

1.3807  X  10-16 

erg/deg  (K) 

Elementary  charge 

e 

4.8032  x  10-10 

statcoulomb 

Electron  mass 

me 

9.1094  x  10-28 

(statcoul) 

g 

Proton  mass 

mp 

1.6726  X  10-24 

g 

Gravitational  constant 

G 

6.6726  X  10-8 

dyne-cm2/g2 

Planck  constant 

h 

6.6261  X  10-27 

erg-sec 

31* 

3- 

to 

4 

1.0546  X  10-27 

erg-sec 

Speed  of  light  in  vacuum 

c 

2.9979  X  1010 

cm/sec 

Proton/electron  mass 

mp/me 

1.8362  x  103 

ratio 

Electron  charge/mass 

e/rrie 

5.2728  X  1017 

statcoul/g 

ratio 

Rydberg  constant 

D  27 r2me4 

R°°  ~  Ch 3 

1.0974  x  105 

cm- 1 

Bohr  radius 

ao  =  h2 /me2 

5.2918  X  10-9 

cm 

Atomic  cross  section 

nao2 

8.7974  x  10-17 

cm2 

Classical  electron  radius 

re  =  e  /me 

2.8179  x  10-13 

cm 

Thomson  cross  section 

(8tt/3  )r2 

6.6525  X  10-25 

cm2 

Compton  wavelength  of 

h/ mec 

2.4263  X  10_1° 

cm 

electron 

h/mec 

3.8616  X  10-11 

cm 

Fine-structure  constant 

a  =  e2 /he 

7.2974  x  10-3 

First  radiation  constant 

a~  1 

ci  =  27 rhe2 

137.04 

3.7418  x  10-5 

erg-cm2/sec 

Second  radiation 

C2  =  he/  k 

1.4388 

cm-deg  (K) 

constant 

Stefan-Boltzmann 

u 

5.6705  X  10“5 

erg/cm2- 

constant 

Wavelength  associated 

^0 

1.2398  X  10-4 

sec-deg4 

cm 

with  1  eV 
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Physical  Quantity 


Symbol 


Value 


Units 


Frequency  associated 

z-'o 

2.4180  X  1014 

Hz 

with  1  eV 

Wave  number  associated 

k0 

8.0655  X  103 

cm- 1 

with  1  eV 

Energy  associated  with 

1.6022  x  10-12 

erg 

1  eV 

Energy  associated  with 

1.9864  x  10-16 

erg 

1  cm-1 

Energy  associated  with 

13.606 

eV 

1  Rydberg 

Energy  associated  with 

8.6174  x  10“5 

eV 

1  deg  Kelvin 

Temperature  associated 

1.1604  x  104 

deg  (K) 

with  1  eV 

Avogadro  number 

Na 

6.0221  x  1023 

mol-1 

Faraday  constant 

F  =  NAe 

2.8925  X  1014 

statcoul/mol 

Gas  constant 

R  =  NAk 

8.3145  X  107 

erg/deg-mol 

Loschmidt’s  number 

n0 

2.6868  X  1019 

_  G 

cm 

(no.  density  at  STP) 

Atomic  mass  unit 

mu 

1.6605  X  10“24 

g 

Standard  temperature 

To 

273.15 

deg  (K) 

Atmospheric  pressure 

po  =  n0kT0 

1.0133  X  106 

dyne/cm2 

Pressure  of  1  mm  Hg 

1.3332  x  103 

dyne/cm2 

(1  torr) 

Molar  volume  at  STP 

Vo  =  RTo/po 

2.2414  x  104 

cm3 

Molar  weight  of  air 

Mair 

28.971 

g 

calorie  (cal) 

4.1868  X  107 

erg 

Gravitational 

9 

980.67 

cm/sec2 

acceleration 
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FORMULA  CONVERSION8 

Here  a  =  102cmm“l,  (3  =  107ergJ_1,  eo  =  8.8542  X  10“  2Fm“L, 
/x  o  =  47 r  x  10“ 7  Hm“1,c=  (eo/io)_1//2  =  2.9979xl08ms“1,  and  h  =  1.0546  X 

_  O  A 

10  J  s.  To  derive  a  dimensionally  correct  SI  formula  from  one  expressed  in 
Gaussian  units,  substitute  for  each  quantity  according  to  Q  =  kQ,  where  k  is 
the  coefficient  in  the  second  column  of  the  table  corresponding  to  Q  (overbars 

denote  variables  expressed  in  Gaussian  units).  Thus,  the  formula  ao  =  h  /me 
for  the  Bohr  radius  becomes  aao  =  (h(3)2 /[(m(3  /  a2)(e2  a(3  / 47reo)],  or  ao  = 
eoh2  / Tzme2 .  To  go  from  SI  to  natural  units  in  which  h  =  c  =  1  (distinguished 
by  a  circumflex),  use  Q  =  k  —  1Q,  where  k  is  the  coefficient  corresponding  to 
Q  in  the  third  column.  Thus  do  =  47reod2/[(md/c)(e2eodc)]  =  4n/rhe2.  (In 
transforming  from  SI  units,  do  not  substitute  for  eo,  /xq,  or  c.) 


Physical  Quantity 

Gaussian  Units  to  SI 

Natural  Units  to  SI 

Capacitance 

a/47reo 

eo  ~  ’ 

Charge 

(oc(3/ 47re0)1/2 

(e0dc)  “ 1/2 

Charge  density 

(/3/47ra5e0)1/2 

(e0dc)“  1/2 

Current 

(a/3/47re0)1/2 

(ix0/hc)1/2 

Current  density 

((3/4na3e0)1/2 

(ix0/hc)1/2 

Electric  field 

(47T/3e0/a3)1/2 

(e0 /hc)1/2 

Electric  potential 

(47T/3e0/a)1/2 

(e0 /hc)1/2 

Electric  conductivity 

(47re0)_1 

eo  1 

Energy 

(3 

(he)-1 

Energy  density 

(3  /  a3 

(he)-1 

Force 

/3/a 

(he)-1 

Frequency 

1 

e-1 

Inductance 

4tt  6o  /  OL 

Ho  1 

Length 

a 

1 

Magnetic  induction 

(4n[3/a3ixo)1/2 

(, Hohc)-1' 2 

Magnetic  intensity 

(47T/x0(3/a3)1/2 

(ix0/hc)1/2 

Mass 

/3/a2 

c/h 

Momentum 

(3  /  a 

h-1 

Power 

(3 

(he2)-1 

Pressure 

(3/a3 

(he)-1 

Resistance 

47reo  /  a 

(eo/h-o)1//2 

Time 

1 

c 

Velocity 

a 

c-1 
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MAXWELL’S  EQUATIONS 


Name  or  Description 

SI 

Gaussian 

Faraday’s  law 

d& 

V  X  E  = - 

dt 

1  OB 

V  x  E  = - 

c  dt 

Ampere’s  law 

dJD 

V  x  H  =  +  J 

dt 

1  OD  4t t 

V  x  H  =  +  J 

c  dt  c 

Poisson  equation 

<1 

D 

■O 

<1 

D 

4 

[Absence  of  magnetic 
monopoles] 

V  •  B  =  0 

V  •  B  =  0 

Lorentz  force  on 

q  (E  +  v  X  B) 

q  (  E  -| - v  X  B  } 

charge  q 

\  c  J 

Constitutive 

D  =  eE 

D  =  eE 

relations 

B=/xH 

B  =  pH 

In  a  plasma,  p.  ~  p. o  =  47 r  X  10  ‘  Hm  1  (Gaussian  units:  p,  ~  1).  The 
permittivity  satisfies  e  ~  eo  =  8.8542  x  10  — 12  F  m  —  1  (Gaussian:  e  ~  1) 
provided  that  all  charge  is  regarded  as  free.  Using  the  drift  approximation 
vi  =  E  X  B / B2  to  calculate  polarization  charge  density  gives  rise  to  a  dielec¬ 
tric  constant  K  =  e/eo  =  1  +  367r  x  109 p/B2  (SI)  =  1  +  47t pc2 / B2  (Gaussian), 
where  p  is  the  mass  density. 

The  electromagnetic  energy  in  volume  V  is  given  by 


W 


dV (H  •  B  +  E  •  D) 


dV (H  •  B  +  E  •  D) 


Poynting’s  theorem  is 


(SI) 


(Gaussian) . 


dW 

dt 


+ 


N  •  dS 


dVJ  •  E, 


where  S  is  the  closed  surface  bounding  V  and  the  Poynting  vector  (energy  flux 
across  S )  is  given  by  N  =  E  x  H  (SI)  or  N  =  cE  x  H/47T  (Gaussian). 
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ELECTRICITY  AND  MAGNETISM 


In  the  following,  e  =  dielectric  permittivity,  /U  =  permeability  of  conduc¬ 
tor,  n'  =  permeability  of  surrounding  medium,  a  =  conductivity,  f  =  <jo/2it  = 
radiation  frequency,  Km  =  fi/fio  and  ne  =  e/eo-  Where  subscripts  are  used, 
‘1’  denotes  a  conducting  medium  and  ‘2’  a  propagating  (lossless  dielectric) 
medium.  All  units  are  SI  unless  otherwise  specified. 


Permittivity  of  free  space 

Permeability  of  free  space 

Resistance  of  free  space 

Capacity  of  parallel  plates  of  area 
A,  separated  by  distance  d 

Capacity  of  concentric  cylinders 
of  length  l,  radii  a,  b 

Capacity  of  concentric  spheres  of 
radii  a,  b 

Self-inductance  of  wire  of  length 
l,  carrying  uniform  current 

Mutual  inductance  of  parallel  wires 
of  length  l,  radius  a,  separated 
by  distance  d 

Inductance  of  circular  loop  of  radius 
b ,  made  of  wire  of  radius  a, 
carrying  uniform  current 

Relaxation  time  in  a  lossy  medium 

Skin  depth  in  a  lossy  medium 

Wave  impedance  in  a  lossy  medium 

Transmission  coefficient  at 
conducting  surface9 
(good  only  for  T  < Cl) 

Field  at  distance  r  from  straight  wire 
carrying  current  I  (amperes) 

Field  at  distance  z  along  axis  from 
circular  loop  of  radius  a 
carrying  current  I 


e0  =  8.8542  x  10-12  Fm-1 

/i o  =  4-7T  x  10-7Hm-1 
=  1.2566  X  10“6Hm_1 

R0  =  (no/e0)1/2  =  376.73  Q 
C  =  eA/d 

C  =  2n el /  \n(b / a) 

C  =  4:TTeab/(b  —  a) 

L  =  ill 

L  =  (p/Z/47r)  [1  +  4  ln(d/a)] 

L  =  b  <^p./  [ln(86/a)  —  2]  +  p./ 4j> 
r  =  e/cr 

S  =  (2/c Uficr)1/2  =  (7T /p.cr)_1//2 
Z  =  [p./(e  +  icr/uj)]1/2 
T  =  4.22  x  10-4(fnrnlne2/(j)112 

Bq=  ill /2nr  tesla 

=  0.21/ r  gauss  (r  in  cm) 

Bz  =  ixa2I/[2{a2  +  z2)3/2] 


21 


ELECTROMAGNETIC  FREQUENCY/ 
WAVELENGTH  BANDS10 


Designation 

Frequency  Range 

Wavelength  Range 

Lower 

Upper 

Lower 

Upper 

ULF* 

30  Hz 

10  Mm 

VF* 

30  Hz 

300  Hz 

1  Mm 

10  Mm 

ELF 

300  Hz 

3  kHz 

100  km 

1  Mm 

VLF 

3  kHz 

30  kHz 

10  km 

100  km 

LF 

30  kHz 

300  kHz 

1  km 

10  km 

MF 

300  kHz 

3  MHz 

100  m 

1  km 

HF 

3  MHz 

30  MHz 

10  m 

100  m 

VHF 

30  MHz 

300  MHz 

1  m 

10  m 

UHF 

300  MHz 

3  GHz 

10  cm 

1  m 

SHFf 

3  GHz 

30  GHz 

1  cm 

10  cm 

S 

2.6 

3.95 

7.6 

11.5 

G 

3.95 

5.85 

5.1 

7.6 

J 

5.3 

8.2 

3.7 

5.7 

H 

7.05 

10.0 

3.0 

4.25 

X 

8.2 

12.4 

2.4 

3.7 

M 

10.0 

15.0 

2.0 

3.0 

P 

12.4 

18.0 

1.67 

2.4 

K 

18.0 

26.5 

1.1 

1.67 

R 

26.5 

40.0 

0.75 

1.1 

EHF 

30  GHz 

300  GHz 

1  mm 

1  cm 

Submillimeter 

300  GHz 

3  THz 

100  p. m 

1  mm 

Infrared 

3  THz 

430  THz 

700  nm 

100  fx m 

Visible 

430  THz 

750  THz 

400  nm 

700  nm 

Ultraviolet 

750  THz 

30  PHz 

10  nm 

400  nm 

X  Ray 

30  PHz 

3  EHz 

100  pm 

10  nm 

Gamma  Ray 

3  EHz 

100  pm 

In  spectroscopy  the  angstrom  is  sometimes  used  (lA  =  10  8  cm  =  0.1  nm). 
*The  boundary  between  ULF  and  VF  (voice  frequencies)  is  variously  defined. 
fThe  SHF  (microwave)  band  is  further  subdivided  approximately  as  shown.11 
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AC  CIRCUITS 


For  a  resistance  R,  inductance  L,  and  capacitance  C  in  series  with 
a  voltage  source  V  =  Vq  exp(zcjt)  (here  i  =  \J  —  1),  the  current  is  given 
by  I  =  dq/dt ,  where  q  satisfies 


L — ^  +  R—  +  — 
dt2  dt  C 


V. 


Solutions  are  q{t )  =  qs  +  Qt,  I(t)  =  Is  +  It,  where  the  steady  state  is 
Is  =  iujqs  =  V / Z  in  terms  of  the  impedance  Z  =  R  +  i(u>L  —  1/ujC)  and 
It  =  dqt/dt.  For  initial  conditions  q( 0)  =  qo  =  go  +  qs,  /( 0)  =  Jo,  the 
transients  can  be  of  three  types,  depending  on  A  =  R2  —  4L/C: 

(a)  Overdamped,  A  >  0 


Io+7+Qo  ,  .v  Io+^-qo  .  . 

qt  =  - exp(— 7_t) - exp(  — 7_|_t) , 

7+  -  7-  7+  -  7- 

7+  (do  +  7- qo)  (  .v  7-  (-^o  +  7+<?o)  ,  .v 

It  =  - exp(  — 7_|_t) - exp(  — 7_t), 

7+  -  7-  7+  -  7- 

where  7-1-  =  (R  ±  A1//2)/2L; 

(b)  Critically  damped,  A  =  0 


gt  =  [go  +  (-To  +  lRQo)t]  exp(— 7 Rt), 

It  =  [Io  -  (/o  +  7J?go)7fi^]  exp(— 7h£), 

where  7#  =  R/2L; 

(c)  Underdamped,  A  <  0 


Qt 

It 


exp(  — 7  Rt), 


'iRQO  +  Io  . 

- sin  <jj\t  +  go  cos  uqt 

vi 

T  ,  (^l2  +  7-R2)go  +  7-R-fo  .  , 

Jo  cos  util  —  - sin(aqt) 

UJl 


exp(-7  Rt), 


Here  u)\  =  ojo (1  —  i?2C'/4L)1/2,  where  cjo  =  (AC)-1/2  is  the  resonant 
frequency.  At  uj  =  luo,  Z  =  R.  The  quality  of  the  circuit  is  Q  =  loqL/ R. 
Instability  results  when  L,  R,  C  are  not  all  of  the  same  sign. 
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DIMENSIONLESS  NUMBERS  OF  FLUID  MECHANICS12 


Name(s) 

Symbol 

Definition 

Significance 

Alfven, 

Karman 

Al,  Ka 

VA/V 

*  (Magnetic  force/ 
inertial  force) 1/7 2 

Bond 

Bd 

(p'  -p)L2g/E 

Gravitational  force/ 
surface  tension 

Boussinesq 

B 

V  /  (2gR)1//2 

(Inertial  force/ 

gravitational  force) 1/72 

Brinkman 

Br 

pV2/kAT 

Viscous  heat /conducted  heat 

Capillary 

Cp 

pV/  E 

Viscous  force/surface  tension 

Carnot 

Ca 

(T2  -Ti)/T2 

Theoretical  Carnot  cycle 
efficiency 

Cauchy, 

Hooke 

Cy,  Hk 

Pv2  /  r  =  m2 

Inertial  force/ 

compressibility  force 

Chandra¬ 

sekhar 

Ch 

B2L2 / pup 

Magnetic  force/dissipative 
forces 

Clausius 

Cl 

LV3  p/kAT 

Kinetic  energy  flow  rate/heat 
conduction  rate 

Cowling 

C 

( VA/V )2  =  Al2 

Magnetic  force/inertial  force 

Crispation 

Cr 

/xk/SL 

Effect  of  diffusion/effect  of 
surface  tension 

Dean 

D 

D3/2V/u(2r)1/2 

Transverse  flow  due  to 

curvature/longitudinal  flow 

[Drag 

coefficient] 

cD 

( P '  ~  p)Lg/ 
pv2 

Drag  force/inertial  force 

Eckert 

E 

V2  /  cpAT 

Kinetic  energy/change  in 
thermal  energy 

Ekman 

Ek 

( u/2QL2)1/ 2  = 
(Ro/Re) 1/2 

(Viscous  force/Coriolis  force) 1/72 

Euler 

Eu 

Ap/pV2 

Pressure  drop  due  to  friction/ 
dynamic  pressure 

Froude 

Fr 

V/igL)1'2 

V/NL 

/(Inertial  force/gravitational  or 
buoyancy  force) 1/72 

Gay-Lussac 

Ga 

1/(3  AT 

Inverse  of  relative  change  in 
volume  during  heating 

Grashof 

Gr 

gL3  (3AT  /  v2 

Buoyancy  force/ viscous  force 

[Hall 

coefficient] 

Ch 

A /rL 

Gyrofrequency  / 

collision  frequency 

*(j)  Also  defined  as  the  inverse  (square)  of  the  quantity  shown. 
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Name(s) 

Symbol 

Definition 

Significance 

Hartmann 

H 

BL/far,)1'2  = 
(Rm  ReC)1/2 

(Magnetic  force/ 

dissipative  force) 1//2 

Knudsen 

Kn 

X/L 

Hydrodynamic  time/ 
collision  time 

Lewis 

Le 

k/V 

*Thermal  conduction/molecular 
diffusion 

Lorentz 

Lo 

V/c 

Magnitude  of  relativistic  effects 

Lundquist 

Lu 

HoLVa/t]  = 

A1  Rm 

J  X  B  force/resistive  magnetic 
diffusion  force 

Mach 

M 

V/Cs 

Magnitude  of  compressibility 
effects 

Magnetic 

Mach 

Mm 

v/vA  =  Ar1 

(Inertial  force/magnetic  force) 1/7 2 

Magnetic 

Reynolds 

Rm 

poLV/rj 

Flow  velocity/magnetic  diffusion 
velocity 

Newton 

Nt 

F/pL2V 2 

Imposed  force/inertial  force 

Nusselt 

N 

aL/k 

Total  heat  transfer/thermal 
conduction 

Peclet 

Pe 

LV/k 

Heat  convection/heat  conduction 

Poisseuille 

Po 

D2  Ap/ pLV 

Pressure  force/viscous  force 

Prandtl 

Pr 

V  /  K 

Momentum  diffusion/ 
heat  diffusion 

Rayleigh 

Ra 

gH3(3AT/vn 

Buoyancy  force/diffusion  force 

Reynolds 

Re 

LV/u 

Inertial  force/ viscous  force 

Richardson 

Ri 

(NH/AV)2 

Buoyancy  effects/ 
vertical  shear  effects 

Rossby 

Ro 

V/2HL  sin  A 

Inertial  force/Coriolis  force 

Schmidt 

Sc 

v/V 

Momentum  diffusion/ 
molecular  diffusion 

Stanton 

St 

a/pcpV 

Thermal  conduction  loss/ 
heat  capacity 

Stefan 

Sf 

crLT3 /k 

Radiated  heat/conducted  heat 

Stokes 

S 

v/L2f 

Viscous  damping  rate/ 
vibration  frequency 

Strouhal 

Sr 

fL/V 

Vibration  speed/flow  velocity 

Taylor 

Ta 

(2  ViL2  lv)2 
R1/2(Ai?)3/2 
■(n/u) 

Centrifugal  force/ viscous  force 
(Centrifugal  force/ 
viscous  force)1//2 

Thring, 

Boltzmann 

Th,  Bo 

pCpV/  ecrT3 

Convective  heat  transport/ 
radiative  heat  transport 

Weber 

W 

pLV2 /T, 

Inertial  force/surface  tension 
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Nomenclature: 

B 

Magnetic  induction 

Ca,  c 

Speeds  of  sound,  light 

Cp 

Specific  heat  at  constant  pressure  (units  ms  K 

D  =  2R 

Pipe  diameter 

F 

Imposed  force 

f 

Vibration  frequency 

9 

Gravitational  acceleration 

H,  L 

Vertical,  horizontal  length  scales 

k  =  pCpK, 

Thermal  conductivity  (units  kgm“  1  s~2) 

N=(g/H)1'2 

Brunt —Vaisala  frequency 

R 

Radius  of  pipe  or  channel 

r 

Radius  of  curvature  of  pipe  or  channel 

VL 

Larmor  radius 

T 

Temperature 

V 

Characteristic  flow  velocity 

VA  =  B/{ix0p)1/ 2 

Alfven  speed 

dT 

Newton’s-law  heat  coefficient,  k -  =  a. AT 

dx 

a. 

(3 

Volumetric  expansion  coefficient,  dV/V  =  (3dT 

r 

Bulk  modulus  (units  kgm~  1  s_  2) 

A R,  AV,  A p,  AT 

Imposed  differences  in  two  radii,  velocities, 
pressures,  or  temperatures 

e 

Surface  emissivity 

r, 

Electrical  resistivity 

K,  T> 

Thermal,  molecular  diffusivities  (units  m2  s  —  1) 

A 

Latitude  of  point  on  earth’s  surface 

A 

Collisional  mean  free  path 

p  =  pu 

Viscosity 

To 

Permeability  of  free  space 

V 

Kinematic  viscosity  (units  m2  s~  L) 

p 

Mass  density  of  fluid  medium 

p' 

Mass  density  of  bubble,  droplet,  or  moving  object 

S 

Surface  tension  (units  kgs-2) 

cr 

Stefan— Boltzmann  constant 

n 

Solid-body  rotational  angular  velocity 
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SHOCKS 


At  a  shock  front  propagating  in  a  magnetized  fluid  at  an  angle  9  with 
respect  to  the  magnetic  induction  B,  the  jump  conditions  are  13,14 

(1)  pU  =  pU  =  q ; 

(2)  pU 2  +  p  +  B j2 /2p  =  pU 2  +  p  +  B  2 /2p,; 

(3)  pUV  -  £> ||  B± / p  =  pUV  -  BllB±/p- 

(4)  B\\  =  % 

(5)  -  VB||  =  UB±  -  VB\\\ 

(6)  ±(U2  +V2)  +  w  +  (UB2  -VB^B^/ppU 

=  i(C/2  +  V2)  +  tD  +  {UB2  -  VB\\Bj_)/ ppU. 

Here  U  and  V  are  components  of  the  fluid  velocity  normal  and  tangential  to 
the  front  in  the  shock  frame;  p  =  1/v  is  the  mass  density;  p  is  the  pressure; 
B±  =  BsinO,  £>||  =  B  cos  9;  p  is  the  magnetic  permeability  (p  =  4-rr  in  cgs 
units);  and  the  specific  enthalpy  is  w  =  e  +  pv ,  where  the  specific  internal 
energy  e  satisfies  de  =  T ds  —  pdv  in  terms  of  the  temperature  T  and  the 
specific  entropy  s.  Quantities  in  the  region  behind  (downstream  from)  the 
front  are  distinguished  by  a  bar.  If  B  =  0,  thenlj 

(7)  U  -U  =  [{p-p){v  -v))1' 2; 

(8)  (p  -  p){v  -  D)-1  =  q2\ 

(9)  w  —  w  =  \  (p  —  p)  (v  +  v) ; 

(10)  e  —  e  =  77  {p  +  p){v  —  v). 

In  what  follows  we  assume  that  the  fluid  is  a  perfect  gas  with  adiabatic  index 
7  =  1  +  2/n,  where  n  is  the  number  of  degrees  of  freedom.  Then  p  =  pRT/m, 
where  R  is  the  universal  gas  constant  and  m  is  the  molar  weight;  the  sound 
speed  is  given  by  Cs 2  =  ( dp/dp)s  =  7 pv,  and  w  =  7e  =  7?w/(7  —  !)•  F°r  a 
general  oblique  shock  in  a  perfect  gas  the  quantity  X  =  r-1  (U /Va)2  satisfies14 

(11)  (A  — /3/a)  (A  — cos2  9)2  =  A  sin2  d<^[l  +  (r—  l)/2a]A  —  cos2  9  where 
r  =  p/p,  a  =  ^[7  +  1  —  (7—  1  )r],  and  (3  =  C  2 /Va2  =  An^p/ B2 . 

The  density  ratio  is  bounded  by 

(12)  1  <  r  <  (7  +  l)/(7  -  1). 

If  the  shock  is  normal  to  B  (i.e.,  if  9  =  7t/2),  then 

(13)  U2  =  (r/a)  {Cs2  +  Aa2  [1  +  (1  -  l/2)(r  -  1)]}; 

(14)  U/U  =  B/B  =  r; 
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(15)  V  =  V; 

(16)  p  =  p  +  (1  —  r~1)pU2  +  (1  —  r2)B2 / 2/i. 

If  0  =  0,  there  are  two  possibilities:  switch-on  shocks,  which  require  (3  <  1  and 
for  which 

(17)  U2  =  rVA2; 

(18)  U  =  V A  2 /U\ 

(19)  B2  =  2B^2{r  —  l)(cr  —  /3); 

(20)  V  =  UBa/B || ; 


(21)  p  =  p  +  pU2{l  -  a.  +  (3){l  -  r  1 ), 

and  acoustic  (hydrodynamic)  shocks,  for  which 

(22)  U 2  =  (r / a)Cs2] 

(23)  U  =  U/r ; 

(24)  V  =  B±  =  0; 

(25)  p  =  p  +  pU2{  1  -  r—  1 ). 

For  acoustic  shocks  the  specific  volume  and  pressure  are  related  by 

(26)  v/v  =  [(7  +  1  )p  +  (7  -  1  )P\  /  [(7  -  1  )p  +  (7  +  1  )P]- 
In  terms  of  the  upstream  Mach  number  M  =  U/Cs , 

(27)  p/p  =  v/v  =  U/U  =  (7  +  l)M2/[(7  -  1)M2  +  2]; 

(28)  p/p  =  (27 M2  -  7  +  1)/ (7  +  1); 

(29)  T/T  =  [(7  -  1  )M2  +  2] (2 7M2  -  7  +  l)/(7  +  1)2M2; 

(30)  M2  =  [(7  -  1)M2  +  2]/[27M2  -7+I]. 

The  entropy  change  across  the  shock  is 

(31)  A s  =  s  -  s  =  cv  In [(p/p)(p/py], 

where  cv  =  R/( 7  —  1  )m  is  the  specific  heat  at  constant  volume;  here  R  is  the 
gas  constant.  In  the  weak-shock  limit  ( M  — >  1), 


(32)  As 


27(7  ~ 

3(7+1) 


(M2  -  1); 


I67  R 
3(7  +  l)m 


(M  —  1) 


The  radius  at  time  t  of  a  strong  spherical  blast  wave  resulting  from  the  explo¬ 
sive  release  of  energy  £  in  a  medium  with  uniform  density  p  is 


(33)  Rs  =  C0(Et2/p )1/5, 


where  Cq  is  a  constant  depending  on  7.  For  7  =  7/5,  Cq  =  1.033. 
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FUNDAMENTAL  PLASMA  PARAMETERS 

All  quantities  are  in  Gaussian  cgs  units  except  temperature  (T,  Te,  T?; ) 
expressed  in  eV  and  ion  mass  (mj)  expressed  in  units  of  the  proton  mass, 
ji  =  mi/mp;  Z  is  charge  state;  k  is  Boltzmann’s  constant;  K  is  wavenumber; 
7  is  the  adiabatic  index;  In  A  is  the  Coulomb  logarithm. 

Frequencies 

electron  gyrofrequency  fce  =  ooce/2Ti  =  2.80  X  106R  Hz 

ujce  =  e£? / mec  =  1.76  X  10  ‘  B  rad/sec 
ion  gyrofrequency  fci  =  Loci/2n  =  1.52  x  103 Z fi~ 1 B  Hz 

u uci  =  ZeB / rriiC  =  9.58  X  10 3 Z/j,~  1 B  rad/sec 
electron  plasma  frequency  fpe  =  cjpe/27r  =  8.98  X  103ne1^2  Hz 

<• ope  =  (47rnee2/me)1/2 

=  5.64  x  10 4ne1//2  rad/sec 
ion  plasma  frequency  fpi  =  upi/ 2n 

=  2.10  X  102Z/i”1/2n*1/2  Hz 

Upi  =  (47miZ2e2/mi)1/2 

=  1.32  x  103Z/i_1/2ni1/2rad/sec 
electron  trapping  rate  VTe  —  (eKE /me)1//2 

=  7.26  X  108K1/2£1/2  sec”1 
ion  trapping  rate  uti  =  (ZeK  E /rrii)1^2 

=  1.69  X  lO7^1/2^:1/2#1/2^-1/2  sec-1 
electron  collision  rate  ve  =  2.91  X  10”  6ne  In  ATe” 3//2  sec- 1 

ion  collision  rate  rq  =4.80  X  10”  8 Z4 /u,~  1//2rq  In  ATi ”  3//2  sec- 1 

Lengths 

electron  deBroglie  length  A  =  h/ (mefcTe)1,/2  =  2.76  X  10”  Te“!//2cm 

classical  distance  of  e2 /kT  =  1.44  x  10’  'T_1  cm 

minimum  approach 

electron  gyroradius  re  =  ^Te/^ce  =  2.38Te  1//2B” 1  cm 

ion  gyroradius  77  =  VTi/^ci 

=  1.02  x  102 ix1/2Z-1Ti1/2B~1  cm 
electron  inertial  length  c/u:pe  =  5.31  X  105ne_1^2cm 

ion  inertial  length  c/ujpi  =  2.28  X  107  (n/ rii)1/2  cm 

Debye  length  Xu  =  (fcT/47me2)1//2  =  7.43  X  102 (T/n)1/2  cm 
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Velocities 

electron  thermal  velocity 
ion  thermal  velocity 
ion  sound  velocity 

Alfven  velocity 

Dimensionless 

(electron/proton  mass  ratio) 1//2 

number  of  particles  in 
Debye  sphere 

Alfven  velocity /speed  of  light 

electron  plasma/gyrofrequency 
ratio 

ion  plasma/gyrofrequency  ratio 
thermal/magnetic  energy  ratio 
magnetic/ion  rest  energy  ratio 

Miscellaneous 

Bohm  diffusion  coefficient 

transverse  Spitzer  resistivity 


vTe  =  ( kTe/me)1// 2 

=  4.19  X  10 7Te1/2  cm/sec 
VTi  =  ( kTi/mi)1/ 2 

=  9.79  X  105n~1/2Ti1/2  cm/sec 

Cs  =  (7 ZkTe/mi )^2 

=  9.79  X  105(7ZTe//u)1/2  cm/sec 
va  =  B  /  (4:irnimi)1/2 

=  2.18  X  101V_1/2^i_1/2Scm/sec 

(■ me/mp)1/ 2  =  2.33  X  10-2  =  1/42.9 
(47r/3)nAD3  =  1.72  x  109T3/2n“1/2 

va/c  =  7.28 ii~  l/2rii~  '  ^ 2 B 
CJpe/cJce  =  3.21  x  10  ~3ne1/2  B~  ' 

Upi/tUci  =  0.137/u1/2ni1/2B_  ! 

(3  =  8tt nkT/B2  =  4.03  X  10-11nTB-2 
B2 /SnnirriiC2  =  26.5  n~1rii~1  B2 

Db  =  {ckT/lQeB) 

=  6.25  X  106TB_1  cm2 /sec 
7U  =  1-15  X  10_14ZlnAT“3/2  sec 
=  1.03  X  10“2Zln  AT“3/2  ficm 


The  anomalous  collision  rate  due  to  low-frequency  ion-sound  turbulence  is 
v*  «  u;peW/kT  =  5.64  x  104ne1/2W/kT  see- 1 , 


where  W  is  the  total  energy  of  waves  with  c u/K  <  VTi- 
Magnetic  pressure  is  given  by 

Pmag  =  B2 / 8n  =  3.98  X  106  (B / Bo)2  dynes/cm2  =  3.93(B / Bo)2  atm, 
where  Bo  =  10  kG  =  1  T. 

Detonation  energy  of  1  kiloton  of  high  explosive  is 

VbkT  =  1012  cal  =  4.2  x  1019  erg. 
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PLASMA  DISPERSION  FUNCTION 


Definition16  (first  form  valid  only  for  Im(  >  0): 

r+oo  eXp  (  —  t2^) 

ry  f  /-\  _  — 1/2  /  V  /  _  o„' _ _ 


Z(  0  =  *• 


t-C 


=  2/  exp  ^  — C2^)  /  dt  exp  ^  — t2j  . 


Physically,  £  =  cc  +  iy  is  the  ratio  of  wave  phase  velocity  to  thermal  velocity. 
Differential  equation: 


—  =-2(l  +  C Z),  Z (0)  =  in1/2; 
d( 


d2Z  dZ 

- h  2C - h  2Z  =  0. 

d(2  d( 


Real  argument  (y  =  0): 


Z (x)  =  exp  ^  —  x2^j  ^in1/2  —  2  J  dt  exp  • 


Imaginary  argument  (x  =  0): 


Z(iy)  =  in1/2  exp  (y2)  [1  -  erf(y)]  . 


Power  series  (small  argument): 


Z{ C)  =  exp  (-C2)  -  2C  (l  -  2C2/3  +  4C4/15  -  8C6/105  +  •■•). 

Asymptotic  series,  |£|  1  (Ref.  IT): 

Z( C)  =  ^7^1/2aexp  (-C2)  -  C_1  (l  +  1/2C2  +  3/4C4  +  15/8C6  +  •••), 

where 

{0  y  >  |  x  |  “  1 

1  \y I  <  M“' 

2  y  <  —  |x|  ~  1 

Symmetry  properties  (the  asterisk  denotes  complex  conjugation): 

Z(C*)  =  -[Z(- C)]*; 

z( c*)  =  [Z(0]  *  +  2/tt1/2  exp[-(C*)2]  (y  >  0). 

Two-pole  approximations18  (good  for  £  in  upper  half  plane  except  when  y  < 
7r 1  ^  2cc2  exp(  — a:2),  a:  A>  1): 

,  x  0.50  +  0.81/  0.50  -  0.81/ 

Z( C)« - ,  a  =  0.51  -  0.81/; 

a-C  a*  +  C 

.  x  0.50  +  0.96/  0.50  -  0.96/ 

Z  ( 0  «  — ; - -x - I - — - — ,  b  =  0.48  —  0.91/. 

fe-c2  fo*+C2 
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COLLISIONS  AND  TRANSPORT 


Temperatures  are  in  eV;  the  corresponding  value  of  Boltzmann’s  constant 
is  k  =  1.60  X  10-12erg/eV;  masses  /i,  p!  are  in  units  of  the  proton  mass; 
ea  =  Zae  is  the  charge  of  species  a.  All  other  units  are  cgs  except  where 
noted. 

Relaxation  Rates 

Rates  are  associated  with  four  relaxation  processes  arising  from  the  in¬ 
teraction  of  test  particles  (labeled  a)  streaming  with  velocity  va  through  a 
background  of  field  particles  (labeled  (3): 


slowing  down 
transverse  diffusion 
parallel  diffusion 
energy  loss 


dva 

dt 

d  i 

Tt{Va 

d  i 

H(v°‘ 

d  o 


a.  |  (3 


_  \  2  a  I  (3  2 

Va)j_  =  V  '  Va 


-  \  2  a  p 

-  Va),|  =  U..  V 


a.  |  /3  _  2 
"'a 


dt 


■Vc.  = 


a  1/3  2 


where  va  =  |vq,|  and  the  averages  are  performed  over  an  ensemble  of  test 
particles  and  a  Maxwellian  field  particle  distribution.  The  exact  formulas  may 
be  written19 


where 


ol  |  /3 

us  IP  = 

=  2 

= 


(1  +  mot/mf3)ilj(xa\(3)v^13; 

(1  —  1/2  xa^)ip(xa^f3)  +  ip'  (xa^) 


v 


a|/3. 

0  ’ 


V 


a|/3. 


uaW  =  2 


ip(xa^)/xa^ 

_l 

(mol/mp)ip(xa^f3)  —  ip'  (xa^) 


v 


ct\(3 
0  ’ 


V 


a|/3 


2  2  a 

=  47rea  ep  X apnp/ma  v 


3. 

ex.  j 


X 


=  mpVa2 /2kTp; 


ip(x) 


dt  t 


1/2  -t. 


ip'  (x) 


dip 
dx  ’ 


and  Xap  =  lnAajg  is  the  Coulomb  logarithm  (see  below).  Limiting  forms  of 
v s ,  u±  and  are  given  in  the  following  table.  All  the  expressions  shown 
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have  units  cm3  sec-  .  Test  particle  energy  e  and  field  particle  temperature  T 
are  both  in  eV;  /i  =  rrii/mp  where  mp  is  the  proton  mass;  Z  is  ion  charge 
state;  in  electron— electron  and  ion— ion  encounters,  field  particle  quantities  are 
distinguished  by  a  prime.  The  two  expressions  given  below  for  each  rate  hold 

for  very  slow  (xa^  <C  1)  and  very  fast  (xa^  /$>  1)  test  particles,  respectively. 


Electron— ion 

r2 


ves'1 /mZ  Xei 
/mZ2 \ei 
u^z/mZ2xei 


Slow 

Fast 

Electron— electron 

el  e  /  x 

/ne  Aee 

«  5.8  x 

10-6T-3//2 

— >  7.7  x 

10-6£-3/2 

el  e  /  x 

/  neXee 

«  5.8  x 

10-6T-1/2e-1 

— >  7.7  x 

10-«e-3/2 

el  e  /  x 

/neXee 

«  2.9  x 

10-6T-1/2e-1 

— ■>  3.9  x 

10~6Te~5/2 

0.23  h3/2T~3/2 

2.5  X  10-V1/2T-1/2e-1 

1.2  x  10-V1/2T-1/2e-1 


->  3.9  X  10-6e-3/2 
->  7.7  x  10-6e-3/2 
->  2.1  x  10-9p.-1Te  — 5/2 


Ion— electron 

j2 


vj  / neZ  Xie 
/neZ2  Xie 

7  I  6  /  ry  2  x 

z^ll  /nez  Xie 

Ion  ion 


v 


i  I  i' 


nvZ*Z'*  A„, 


1.6  x  10-V-1T-3/2 
3.2  x  10-9p.-1T-1/2e-1 
1.6  X  10-9/i-1T-1/2e-1 


6.8  X  10 


-8  V 


n  /  2 


-  1.7  X  10-4m1/26-3/2 
,i.8xio-yl/2r3/2 
-+  1.7  X  10_<V/2Te_5/2 


U 

1  +  ~ 
I1  . 


-1/2 


t—  3/2 


->  9.0  x  10 


-8 


- 1 - 7 

/i  /r 


/i 


1/2 


=  3/2 


V 


nt’Z2Z'l\ti, 


i  A  xx  in~7  'l/2  — 1/Ti— 1/2  -1 

1.4  x  10  ii  '  T  '  e 


n.8xio-yl/2r3/2 


I/,, 


™,'Z2Z'2Ai4, 


xi  q  . ,  in  —  ^  / 1/2  —  1  rji  —  1/2  —  1 

6.8  X  10  /i  '  /x  T  '  e 


n  n  xx  in-8  l/2  '-1^  -5/2 

— >  9.0  X  10  ii  1  ji  Te  ' 

In  the  same  limits,  the  energy  transfer  rate  follows  from  the  identity 

ve  =  2vs  -  -  v\\ , 

except  for  the  case  of  fast  electrons  or  fast  ions  scattered  by  ions,  where  the 
leading  terms  cancel.  Then  the  appropriate  forms  are 


zy 


->  4.2  x  10  9 rin  Z2 Xf 


—  3/2  —  i 


e  -'-fj,  -  -  8.9  X  10  (/i/T)~'  e 


1/2-1 


exp(—  1836p.e/T) 


sec 


-l 
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and 


ui]i'  — >  1.8  X  10~7 ni/  Z2 Z'2 Xii/ 


-3/2  1/2 


e  '  /x  '  / fx  —  1.1  (///T)  '  e  exp(  —  /x'e/T) 


sec 


-l 


In  general,  the  energy  transfer  rate  is  positive  for  e  >  6^*  and  nega¬ 

tive  for  e  <  ea*,  where  x*  =  (mp /ma)ea* /T^ 3  is  the  solution  of  t//(x*)  = 
(mo,  |mjg)'0(x*).  The  ratio  eQ  */Tp  is  given  for  a  number  of  specific  a,  f3  in  the 
following  table: 


CK  |/3 

i\ 

e  e 

e,  t 

i  e|p  e 

D  e 

T,  e 

He3  e 

He4 

Tfi 

1.5  1 

0.98 

4.8  x  10“3  2.6  x 

10-3  1.8  x  10-3  1.4  x  10-3 

When  both  species  are  near  Maxwellian,  with  Tj  A  Te,  there  are  just 
two  characteristic  collision  rates.  For  Z  =  1, 


r'e  =  2.9  X  10  6n\Te  3 ^  sec  1; 
Vi  =  4.8  X  10_8nAT^  _3/2p._  1/2  s 


Temperature  Isotropization 

Isotropization  is  described  by 


dT±  _  1  dT || 

dt  2  dt 


-i4(Tl  -  T||), 


where,  if  A  =  T±/T^  —  1  >  0, 


OL 


V rj-1 


2  2 

2y/ne0,  ep  TictXap  ^  —  2 

ma1/2(fcT|1)3/2 


-3+  (A +  3) 


tan  1(t41//2) 
A1/2 


If  A  <  0,  tan  1  (A1//2)/A1//2  is  replaced  by  tanh  (  —  A)1/2/(  —  A)1/2.  For 
Tj_  «  T||  =  T, 

^  =  8.2  x  10  7 n\T  3/72  sec  1; 

vlp  =  1.9  X  10  8nXZ2p  1/2t  3/72  sec  1. 
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Thermal  Equilibration 

If  the  components  of  a  plasma  have  different  temperatures,  but  no  rela¬ 
tive  drift,  equilibration  is  described  by 


dTa 

dt 


^Vl/3(Tg  -T«), 
0 


where 

pcd/3  _  1  g  x  10-19 
For  electrons  and  ions  with  Te  : 


(mo,  mp ) 1//2  2  Zp  2np  Xap 

{maTp  +  mpTa)3/2 

a  Ti  =  T ,  this  implies 


-l 

sec 


u^z /rii  =  /ne  =  3.2  x  10  9Z2 A/ pT3 7/2 cm3  sec  1. 


Coulomb  Logarithm 

For  test  particles  of  mass  mo,  and  charge  ea  =  scattering  off  field 

particles  of  mass  mp  and  charge  ep  =  Zpe,  the  Coulomb  logarithm  is  defined 

as  A  =  In  A  =  ln(rmax/rmin).  Here  rmin  is  the  larger  of  eaep /mapu2  and 
h/2molpu,  averaged  over  both  particle  velocity  distributions,  where  map  = 
mamp/(ma  +mp)  and  u  =  va  -  v^;  rmax  =  (47 r  ^  n7e72//cT7)_1/2,  where 
the  summation  extends  over  all  species  7  for  which  u2  <  vt- y2  =  fcT7/m7.  If 
this  inequality  cannot  be  satisfied,  or  if  either  'UCJcq,-1  <  rmax  or  ,uc<7Cjg_1  < 
rmax)  the  theory  breaks  down.  Typically  A  ~  10—20.  Corrections  to  the  trans¬ 
port  coefficients  are  0( A~  L);  hence  the  theory  is  good  only  to  ~  10%  and  fails 
when  A  \  # 

The  following  cases  are  of  particular  interest: 

(a)  Thermal  electron— electron  collisions 

Aee  =  23  -  In (ne1/2Te-3/2),  Te  <  lOeV; 

=  24  -  ln^e172^-1),  Te  >  lOeV. 


(b)  Electron— ion  collisions 

Xei  =  Xie  =23 -In  [ne1/2ZT~3/2y 
=  24  -  In  (ne172^1), 

=  30  -  In  (rn1/2Ti~3/2 Z2 


Time / rrii  <  Te  <  10 Z2  eV; 
Time / mi  <  10Z2  eV  <  Te 
Te  <  TiZme/mi. 
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(c)  Mixed  ion— ion  collisions 


Xu'  =  Xi'i  =  23  -  In 


ZZ'(n  +  fj,')  / rijZ2  (  nifZ' 

flT-/  +  n'Ti  y  Ti 


2  X  1/2' 


+ 


(d)  Counterstreaming  ions  (relative  velocity  V£>  =  /3_dc)  in  the  presence  of 
warm  electrons,  kTi/rrii,  kTi/ /m^  <  vjj2  <  kTe/me 


Xu'  =  Xi'i  =  35  -  In 


ZZ'  (fj,  +  n')  f  nt 


1/2 


HH'/3d2 


Fokker-Planck  Equation 

Dr  df 


Dt 


dt 


+  v  ■  V/  +  F  •  Vv/  = 


or 

dt 


coll 


where  F  is  an  external  force  field.  The  general  form  of  the  collision  integral  is 

(dr/m)  coU  =  -  v  Vv  ■ ja|/3’ with 


2  2 

t«I/3  n  \  e«  eP  I  1 3  //  2/  \ 

J  1  =  2TcXaf3 -  a  v  (u  I  —  uu)« 

mn 


-3 


(  —  /“(v)Vv,/9(v')  -  —  /'5(v')Vv/“(v)l 

L  m/3  J 


(Landau  form)  where  u  =  v;  —  v  and  /  is  the  unit  dyad,  or  alternatively, 


2  _  2 


Ja|/3  =  l~A,v.j  — — |/“(v)VvH(v)  -  ivv  ■  |/“(v)VvVvG(v) 


}■ 


where  the  Rosenbluth  potentials  are 


G(v)  =  /  /^(v/)wd3r;/ 


if(v)  =  1  + 


m. 


mp 


r  (3  /  /  \  —  1  73  / 

j  [y  )u  a  v  . 
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If  species  a  is  a  weak  beam  (number  and  energy  density  small  compared  with 
background)  streaming  through  a  Maxwellian  plasma,  then 


Ta|/3  Triot  ot\(3  r  ol  1  ck|/3  j-a 

J  = - ; - i's  v/  -  -v  vv  ■  Vv/ 

+  mg  2  II 

ct  |  (3  /  2  /  \  t— r  r  ct 

-  ~  VVJ  ■  Vv/  ■ 

B-G-K  Collision  Operator 

For  distribution  functions  with  no  large  gradients  in  velocity  space,  the 
Fokker-Planck  collision  terms  can  be  approximated  according  to 

Dfe 

=  Vee(Fe  -  fe)  +  Uei(Fe  -  /e); 

Dfi 

=  Vie(Fi  -  fi )  +  Vii{Fi  -  fi). 

The  respective  slowing-down  rates  given  in  the  Relaxation  Rate  section 

above  can  be  used  for  uap,  assuming  slow  ions  and  fast  electrons,  with  e  re¬ 
placed  by  Ta.  (For  uee  and  va,  one  can  equally  well  use  u±,  and  the  result 
is  insensitive  to  whether  the  slow-  or  fast-test-particle  limit  is  employed.)  The 
Maxwellians  Fa  and  Fa  are  given  by 


where  n^,  vQ  and  Ta  are  the  number  density,  mean  drift  velocity,  and  effective 
temperature  obtained  by  taking  moments  of  fa.  Some  latitude  in  the  definition 
of  Ta  and  va  is  possible;20  one  choice  is  Te  =  Ti,  Ti  =  Te,  ve  =  v.j,  Vj  =  ve. 

Transport  Coefficients 

Transport  equations  for  a  multispecies  plasma: 


dana 

dt 


+  na  V 


va 


0; 


d“va 

mana  — - —  =  —  Vpa  —  V  •  Fa  +  Zaent 
dt 


E  +  —  va  x  B 

c 


+  ha ; 
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(J  U/  r\j  -L  q/  _ 

2  n"  — dt - ^  Pa  V  '  Va  =  _V  ’  q“  _  :  Vv“  +  Q®' 

Here  da / dt  =  (3/c3t  +  va  ■  V;  =  nakTa ,  where  /c  is  Boltzmann’s  constant; 
Ra  =  T;  Ra/j  and  Qa/ 3,  where  Ra/3  and  Qa/3  are  respectively 

the  momentum  and  energy  gained  by  the  ath  species  through  collisions  with 
the  /3th;  Pa  is  the  stress  tensor;  and  qa  is  the  heat  flow. 

The  transport  coefficients  in  a  simple  two-component  plasma  (electrons 
and  singly  charged  ions)  are  tabulated  below.  Here  ||  and  _L  refer  to  the  di¬ 
rection  of  the  magnetic  field  B  =  b B\  u  =  ve  —  Vj  is  the  relative  streaming 
velocity;  ne  =  ni  =  n;  j  =  —  neu  is  the  current;  c uce  =  1.76  X  10  B  sec  and 
uJci  =  (nie/mi)u)ce  are  the  electron  and  ion  gyrofrequencies,  respectively;  and 
the  basic  collisional  times  are  taken  to  be 


Te 


3y ^Z{kTef/2 
4\Z2tt  nXe4 


rj~i  3/2 

3.44  x  105  — - - sec, 

n  X 


where  A  is  the  Coulomb  logarithm,  and 


Ti 


3^/Td-(kTi)3/2 

AQnn  Ae4 


2.09  x  107 


rji  3/2 
-L  i 

nX 


sec. 


In  the  limit  of  large  fields  (uiCaTa  /$>  1 ,  a  =  i,  e)  the  transport  processes  may 

2 1 

be  summarized  as  follows: 


momentum  transfer 
frictional  force 

electrical 

conductivities 

thermal  force 
ion  heating 


R-ei —  Rie  —  R  —  Ru  "f  R-T1 1 
Ru  =  ne(  j||/<T||  +jx/crx); 

2 

cr ||  =  1.96cr_|_;  cr±  =  ne  re/me; 


3n 

Rt  =  —  0.7 In V i|  ( kTe )  —  - b  x  V_l  ( kTe)\ 

2c L>ceTe 

3 me  nk 

Qi  = - (' Te-Ti ); 

mi  re 


electron  heating 
ion  heat  flux 

ion  thermal 
conductivities 
electron  heat  flux 

frictional  heat  flux 


Qe 

q  i 

i 

K\\ 

qe 


-Qi  -  R  u; 


—  «||  V ||  (kTQ  -  V i  (kTQ  +  n\h  X  V_L(fcTi); 
nkTiTi  •  2  nkTi  •  5  nkTi 

Q  O  .  Kl  —  _ _.  Kl  —  _ 

O.J  ,  n_,_  —  ,  KA  —  , 

mi  miLujTi  2  rniUJd 

e  ,  e 

qu  +  q  t; 

Q/n  U'  I  1 

0.71n/cTeU||  H - —  b  x  uj_; 

2cuce  Te 
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thermal  gradient 
heat  flux 


=  —  rcjjV||  (kTe)  -  Ke±V ±(kTe)  —  K®b  x  Vj _(kTe); 


nkTpTP 


electron  thermal  /«jj  =  3.2 - ;  /-c^  =  4.7 

conductivities  me 

stress  te 
species) 


nkTp 


meuJc^Te  ’ 


e 

ka  = 


5  nkTe 
2meujc 


r/  o  771 

stress  tensor  (either  Pxx  —  —  —  (Wxx  +  Wyy) - —  (Wxx  —  Wyy )  —  rj3Wxy; 

Pyy  =  - y  (W**  +  WB„)  +  -  Wyy)  +  V3 Wxy-, 

Tj  3 

Pxy —  fyx  —  k!\Wxy  “I"  “T- (l^xsc  IFyy); 

Pxz —  Pzx  —  ?72  Tl4pVyz‘-i 

Pyz —  -Pzy  —  772 hf^yz  4“  TlA^^xz] 

Pzz  =  -rjoWzz 

(here  the  z  axis  is  defined  parallel  to  B); 

3  nkTi 


ion  viscosity 

X 

Vo 

=  0.96  nkTiTi] 

X 

nkTi  i 

v3 

=  ;  Va  = 

CX 

electron  viscosity 

e 

Vo 

=  0.73  n/cTere; 

e 

nkTp  e 

v3 

=  o  ;  va 

^ce 

IOcj  2T;  ’ 


7?2  = 


6  nkTi 
5c<J  2  Ti  ’ 


nkTi 


OJr 


nkTp 


nkTe 

“clTe 


;  vl  =  2-° 


nkTp 

9  ’ 

WceTe 


For  both  species  the  rate-of-strain  tensor  is  defined  as 


Wjk  = 


dvi  dv 

+ 


k 


dxk  dxj  3 


-  -5jk\7  •  v. 


When  B  =  0  the  following  simplifications  occur: 


Ru  =  Txej/ a || ;  RT  =  -0.71nV(fcTe);  q^  =  —  V(fcTj); 

q|j  =  0.71nfcTeu;  q^  =  —  /-cjj  V(kTe);  Pjk  =  —r]0Wjk. 

For  uj pp t p  1  u >dTi,  the  electrons  obey  the  high-field  expressions  and  the 

ions  obey  the  zero-field  expressions. 

Collisional  transport  theory  is  applicable  when  (1)  macroscopic  time  rates 
of  change  satisfy  d/dt  <C  1/r,  where  r  is  the  longest  collisional  time  scale,  and 
(in  the  absence  of  a  magnetic  field)  (2)  macroscopic  length  scales  L  satisfy  L  ^ > 
Z,  where  l  =  vr  is  the  mean  free  path.  In  a  strong  field,  cjcer  1,  condition 
(2)  is  replaced  by  L  y  l  and  L±  \/lre  ( L±  re  in  a  uniform  field), 
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where  L  ||  is  a  macroscopic  scale  parallel  to  the  field  B  and  Lj_  is  the  smaller 
of  B/\X/±B\  and  the  transverse  plasma  dimension.  In  addition,  the  standard 
transport  coefficients  are  valid  only  when  (3)  the  Coulomb  logarithm  satisfies 
A  1;  (4)  the  electron  gyroradius  satisfies  re  Xjj,  or  87 rnemec2  ^$>  B2 ;  (5) 
relative  drifts  u  =  va  —  vp  between  two  species  are  small  compared  with  the 

thermal  velocities,  i.e.,  u 2  <C  kTa/ma,  kTp/mp ;  and  (6)  anomalous  transport 
processes  owing  to  microinstabilities  are  negligible. 

Weakly  Ionized  Plasmas 

Collision  frequency  for  scattering  of  charged  particles  of  species  a.  by 
neutrals  is 

Vol  =  ri0(Js  (feTa/ma)  , 

where  no  is  the  neutral  density  and  cr"\°  is  the  cross  section,  typically  ~ 

_ ip  0 

5  X  10  cm  and  weakly  dependent  on  temperature. 

When  the  system  is  small  compared  with  a  Debye  length,  L  <C  Aq,  the 
charged  particle  diffusion  coefficients  are 

Da  =  kTa/maua, 

In  the  opposite  limit,  both  species  diffuse  at  the  ambipolar  rate 

D  _  k'i^e  —  k'eDi  ( Ti  +  Te)DiDe 

A  m  —  He  TiDe  +  TeDi 

where  Ha  =  &a  /  is  the  mobility.  The  conductivity  <ja  satisfies  era  = 

Tla&a  Ha  • 

In  the  presence  of  a  magnetic  field  B  the  scalars  h  and  a  become  tensors, 
J“  =  <ra  E  =  <r^E||  +  ff“Ei  +  a“E  x  b, 


where  b  =  B  /B  and 


<2”  1 1  —  Tla&a  /'^a^'a'i 

=  0T||  /{ua  +  cvca); 
crA  =  cr i|  Ua^ca/{^a  +^ca)- 


Here  cr j_  and  crA  are  the  Pedersen  and  Hall  conductivities,  respectively. 
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IONOSPHERIC  PARAMETERS 


23 


The  following  tables  give  average  nighttime  values.  Where  two  numbers 
are  entered,  the  first  refers  to  the  lower  and  the  second  to  the  upper  portion 
of  the  layer. 


Quantity 

E  Region 

F  Region 

Altitude  (km) 

90-160 

160-500 

Number  density  (m  —  3) 

1.5  X  1010— 3.0  X  1010 

5  X  1010— 2  x  1011 

Height-integrated  number 

9  X  1014 

4.5  X  1015 

density  (m  —  2) 

Ion-neutral  collision 

2  x  103-102 

0.5-0.05 

frequency  (sec-1) 

Ion  gyro-/collision 

0.09-2.0 

4.6  x  102— 5.0  x  103 

frequency  ratio 

Ion  Pederson  factor 

0.09-0.5 

2.2  x  10-3— 2  x  10-4 

Ki/{  1  “b  Ki  2 ) 

Ion  Hall  factor 

8  X  10-4— 0.8 

1.0 

Kq2/(1  +  /-c^2) 

Electron-neutral  collision 

1.5  x  104— 9.0  x  102 

80-10 

frequency 

Electron  gyro-/collision 

4.1  x  102— 6.9  x  103 

7.8  x  104— 6.2  x  105 

frequency  ratio  ree 

Electron  Pedersen  factor 

2.7  x  10-3— 1.5  x  10-4 

10-5— 1.5  x  10-6 

Ke /(I  +  ^e2) 

Electron  Hall  factor 

1.0 

1.0 

Ke2/ (1  +  Ke2) 

Mean  molecular  weight 

28-26 

22-16 

Ion  gyrofrequency  (sec-1) 

180-190 

230-300 

Neutral  diffusion 

30-5  X  103 

105 

coefficient  (m2  sec- 1) 

The  terrestrial  magnetic  field  in  the  lower  ionosphere  at  equatorial  lati¬ 
tudes  is  approximately  Bq  =  0.35  Xl0  —  4  tesla.  The  earth’s  radius  is  Re  =  6371 
km. 
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SOLAR  PHYSICS  PARAMETERS24 


Parameter 

Symbol 

Value 

Units 

Total  mass 

M0 

1.99  X  1033 

g 

Radius 

Rq 

6.96  X  1010 

cm 

Surface  gravity 

90 

2.74  x  104 

cm  s“" 

Escape  speed 

V  oo 

6.18  x  107 

cm  s—  1 

Upward  mass  flux  in  spicules 

— 

1.6  x  10-9 

g cm-2  s“  ! 

Vertically  integrated  atmospheric  density 

— 

4.28 

g  cm-2 

Sunspot  magnetic  field  strength 

■^max 

2500-3500 

G 

Surface  effective  temperature 

T0 

5770 

K 

Radiant  power 

Cq 

3.83  X  1033 

erg  s“  1 

Radiant  flux  density 

T 

6.28  X  1010 

erg  cm~2s_1 

Optical  depth  at  500  nm,  measured 

T5 

0.99 

— 

from  photosphere 

Astronomical  unit  (radius  of  earth’s  orbit) 

AU 

1.50  X  1013 

cm 

Solar  constant  (intensity  at  1AU) 

/ 

1.36  X  106 

erg  cm-2  s  — 1 

Chromosphere  and  Corona20 


Parameter  (Units) 

Quiet 

Sun 

Coronal 

Hole 

Active 

Region 

Chromospheric  radiation  losses 
(erg cm-2  s-1) 

Low  chromosphere 

2  x  106 

2  x  106 

>  io7 

Middle  chromosphere 

2  x  106 

2  x  106 

107 

Upper  chromosphere 

3  X  105 

3  X  105 

2  x  106 

Total 

4  x  106 

4  x  106 

>  2  x  107 

Transition  layer  pressure  (dyne  cm-2) 

0.2 

0.07 

2 

Coronal  temperature  (K)  at  1.1  Rq 

1. 1-1.6  x  106 

106 

2.5  x  106 

Coronal  energy  losses  (ergcm“  2s_1) 

Conduction 

2  x  105 

6  x  104 

105-107 

Radiation 

105 

104 

5  x  106 

Solar  Wind 

<  5  x  104 

7  x  105 

<  105 

Total 

3  x  105 

8  x  105 

107 

Solar  wind  mass  loss  (gcm_2s_1) 

<  2  x  10-11 

2  x  10_1° 

<  4  x  10-11 
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THERMONUCLEAR  FUSION 


26 


Natural  abundance  of  isotopes: 

hydrogen 

helium 

lithium 


Mass  ratios: 


n D  / n H  =  1-5  X  10  4 
nHe3/™He4  =  1-3  X  10“6 
nLi6/^Li7  =  0-08 


me/mD 


=  2.72  x  10-4  =  1/3670 


(■ me/mD)1/2=  1.65  X  10-2  =  1/60.6 


me/mT  =  1.82  x  10-4  =  1/5496 


(me/mT)1/2=  1-35  X  10-2  =  1/74.1 

Absorbed  radiation  dose  is  measured  in  rads:  1  rad  =  102ergg  —  1.  The  curie 
(abbreviated  Ci)  is  a  measure  of  radioactivity:  1  curie  =  3.7x  1010  counts  sec-  1  . 


Fusion  reactions  (branching  ratios  are  correct  for  energies  near  the  cross  section 
peaks;  a  negative  yield  means  the  reaction  is  endothermic):27 


(la) 

D  +  D  ■ 

5Q^  >T(1.°1  MeV)  +  p(3.02  MeV) 

(lb) 

5Q^  >He3 (0.82  MeV)  +  n(2.45  MeV) 

(2) 

D  +  T 

- >He4(3.5MeV)  +  n(14.1  MeV) 

(3) 

D  +  He3- 

- >He4(3.6MeV)  +  p(14.7MeV) 

(4) 

T  +  T 

- >He4  +  2n  +  11.3  MeV 

(5a) 

He3  +  T 

5lV  >He  "F  P  +  n  +  12.1  MeV 

(5b) 

43^  >He4(4.8  MeV)  +  D(9.5MeV) 

(5c) 

Q0/  >He5 (2.4  MeV)  +  p(11.9MeV) 

(6) 

P  +  Li6  - 

- >He4(1.7MeV)  +  He3(2.3MeV) 

(7a) 

p  +  Li7  - 

-2^2  He4  +  17.3  MeV 

(7b) 

>Be7  +  n  -  1.6  MeV 

(8) 

D  +  Li6  - 

- >2He4  +  22.4  MeV 

(9) 

P  +  B11  - 

- >3  He4  +  8.7  MeV 

(10) 

n  +  Li6  - 

- >He4(2.1  MeV)  +  T(2.7MeV) 

The  total  cross  section  in  barns  (1  barn  =  10  —  24  cm2)  as  a  function  of  E,  the 
energy  in  keV  of  the  incident  particle  [the  first  ion  on  the  left  side  of  Eqs. 
(1)— (5)],  assuming  the  target  ion  at  rest,  can  be  fitted  by28 


cr 7^  (E) 


Or 

+ 

\a4  -  A3E)2  +  l" 

-l 

A2 

E 

exp(AiE  1/2)  —  1 
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where  the  Duane  coefficients  Aj  for  the  principle  fusion  reactions  are  as  follows: 


1 

D-D 

(la) 

D-D 

(lb) 

D-T 

(2) 

D  He3 

(3) 

T-T 

(4) 

T -He3 
(5a— c) 

Ai 

46.097 

47.88 

45.95 

89.27 

38.39 

123.1 

A2 

372 

482 

50200 

25900 

448 

11250 

A3 

4.36  X  10-4 

3.08  X  10-4 

1.368  X  10-2 

3.98  X  10-3 

1.02  x  10-3 

0 

A4 

1.220 

1.177 

1.076 

1.297 

2.09 

0 

A5 

0 

0 

409 

647 

0 

0 

Reaction  rates  crv  (in  cm3  sec  '),  averaged  over  Maxwellian  distributions: 


Temperature 

(keV) 

D-D 

(la  +  lb) 

D-T 

(2) 

D  He3 

(3) 

T-T 

(4) 

T  He3 

(5a— c) 

1.0 

1.5  X  10-22 

5.5  x  10-21 

IQ-26 

3.3  x  10-22 

O 

1 

to 

00 

2.0 

5.4  x  10“21 

2.6  x  10-19 

1.4  x  10-23 

7.1  x  10“21 

IQ-25 

5.0 

1.8  x  10-19 

1.3  x  10-17 

6.7  x  10“21 

1.4  x  10-19 

2.1  x  10-22 

10.0 

1.2  x  10-18 

1.1  x  10-16 

2.3  x  10-19 

7.2  x  10-19 

1.2  x  10_2° 

20.0 

5.2  x  10-18 

4.2  x  10-16 

3.8  x  10-18 

2.5  x  10-18 

2.6  x  10-19 

50.0 

2.1  x  10-17 

8.7  x  10-16 

5.4  x  10-17 

8.7  x  10“18 

5.3  x  10“18 

100.0 

4.5  x  10~17 

8.5  x  10-16 

1.6  x  10-16 

1.9  x  10-17 

2.7  x  10-17 

200.0 

8.8  x  10”17 

6.3  x  10-16 

2.4  x  10-16 

4.2  x  10-17 

9.2  x  10-17 

500.0 

1.8  x  10-16 

3.7  x  10-16 

2.3  x  10-16 

8.4  x  10-17 

2.9  x  10-16 

1000.0 

2.2  x  10-16 

2.7  x  10-16 

1.8  x  10-16 

8.0  x  IQ"17 

5.2  x  10-16 

For  low  energies  (T  A  25  keV)  the  data  may  be  represented  by 

(crv) d d  =  2.33  X  10  14T  2^3  exp(— -18.76T  ^3)cm3  sec  1; 


(orv) £>t  =  3.68  X  io-12T-2/3  exp(  — 19.94T_1//3)  cm3  sec-1 , 
where  T  is  measured  in  keV. 

The  power  density  released  in  the  form  of  charged  particles  is 

Pdd  =  3.3  X  10 ~13 nr)2 (crv) dd  watt  cm-3  (including  the  subsequent 

D— T  reaction); 

Pdt  =  5.6  X  10“  3ri]jnT(crv) dt  watt  cm"  '; 

^DHe^  —  2.9  X  10  watt  cm 
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RELATIVISTIC  ELECTRON  BEAMS 

Here  7  =  (1  —  /32)- 1//2  is  the  relativistic  scaling  factor;  quantities  in 
analytic  formulas  are  expressed  in  SI  or  cgs  units,  as  indicated;  in  numerical 
formulas,  I  is  in  amperes  (A),  B  is  in  gauss  (G),  electron  linear  density  N  is 
in  cm-1,  and  temperature,  voltage  and  energy  are  in  MeV;  (3Z  =  vz/c\  k  is 
Boltzmann’s  constant. 

Relativistic  electron  gyroradius: 

re  =  - (y2  —  1)1//2  (cgs)  =  1.70  x  103(72  —  l)1^2  B~x  cm. 

eB 

Relativistic  electron  energy: 

W  =  me2 7  =  O.5II7  MeV. 

Bennett  pinch  condition: 

1 2  =  2 Nk{Te  +  Ti)c 2  (cgs)  =  3.20  X  10“4A(Te  +  Tz)  A2. 
Alfven-Lawson  limit: 

I A  =  (mc3/e)pz  7  (cgs)  =  (47rmc/p.0e)/3z7  (SI)  =  1.70  X  104/3Z7  A. 

The  ratio  of  net  current  to  I  a  is 

1  _  v 
I A  7 

Here  u  =  Nre  is  the  Budker  number,  where  re  =  e2/mc2  =  2.82  X  10_  13  cm 
is  the  classical  electron  radius.  Beam  electron  number  density  is 

rib  =  2.08  X  108 J (3  1  cm  3, 

where  J  is  the  current  density  in  Acm-2.  For  a  uniform  beam  of  radius  a  (in 
cm), 

rib  =  6.63  X  10 ‘la  2 (3  1  cm  3, 

and 

2 re  _  v 
a  7 
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Child’s  law:  (non-relativistic)  space-charge-limited  current  density  between 
parallel  plates  with  voltage  drop  V  (in  MV)  and  separation  cL  (in  cm)  is 

J  =  2.34  x  103V3/2d-2  A  cm-2. 


The  saturated  parapotential  current  (magnetically  self-limited  flow  along  equi- 
potentials  in  pinched  diodes  and  transmission  lines)  is29 


Ip  =  8.5  X  10" 


Gy  In  7  +  (7  —  1) 


1/2 


A, 


where  G  is  a  geometrical  factor  depending  on  the  diode  structure: 

w  for  parallel  plane  cathode  and  anode 

of  width  w,  separation  d; 

-1 

for  cylinders  of  radii  Ri  (inner)  and  R2  (outer); 

for  conical  cathode  of  radius  Rc,  maximum 
do  separation  do  (at  r  =  Rc )  from  plane  anode. 

For  f3  — >  0  (7  — >  1),  both  I  a  and  Ip  vanish. 

The  condition  for  a  longitudinal  magnetic  field  Bz  to  suppress  filamentation 
in  a  beam  of  current  density  J  (in  Acm~2)  is 

Bz  >  47  (3Z  (7  J)1/2  G. 

Voltage  registered  by  Rogowski  coil  of  minor  cross-sectional  area  A,  n  turns, 
major  radius  a,  inductance  L,  external  resistance  R  and  capacitance  G  (all  in 
SI): 


vr  — 


2nd 


G  =  (in  — 

V  Ri 

/ — 1  _  Rc 


externally  integrated  V  =  (1/ RC)(nAfj,oI /2ira)] 

self-integrating  V  =  (R/ L)(nAnoI /2na)  =  Rl/n. 

X-ray  production,  target  with  average  atomic  number  Z  ( V  <  5  MeV): 
r]  =  x-ray  power/beam  power  =  7  X  10  4 ZV. 

X-ray  dose  at  1  meter  generated  by  an  e-beam  depositing  total  charge  Q 
coulombs  while  V  >  0.84Vmax  in  material  with  charge  state  Z: 

D  =  150Vm2x8QZ1/2  rads. 
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BEAM  INSTABILITIES30 


Name 

Conditions 

Saturation  Mechanism 

Electron- 

electron 

Vd  >  Vej,  j  =  1,2 

Electron  trapping  until 

Cei  -  Vd 

Buneman 

Vd  >  (M/m)1/3W, 

Vd  >  Ve 

Electron  trapping  until 

Ve  -  Vd 

Beam-plasma 

vb  >  (np/nb)V5Vb 

Trapping  of  beam  electrons 

Weak  beam- 
plasma 

Vb  <  (np/nb)1/3V b 

Quasilinear  or  nonlinear 
(mode  coupling) 

Beam-plasma 

(hot-electron) 

Ve>Vb>  Vb 

Quasilinear  or  nonlinear 

Ion  acoustic 

Te  >  Ti  ,  Vd  >  Cs 

Quasilinear,  ion  tail  form¬ 
ation,  nonlinear  scattering, 
or  resonance  broadening. 

Anisotropic 

temperature 

(hydro) 

Te _L  >  2Te|| 

Isotropization 

Ion  cyclotron 

Vd  >  20 Vi  (for 

Te  «  Ti) 

Ion  heating 

Beam-cyclotron 

(hydro) 

vd  >  Cs 

Resonance  broadening 

Modified  two- 
stream  (hydro) 

Vd  <  (1  +  (3)1/2Va, 

vd  >  cs 

Trapping 

Ion-ion  (equal 
beams) 

U  <  2(1  +(3)1/2Va 

Ion  trapping 

Ion-ion  (equal 
beams) 

U  <  2  Cs 

Ion  trapping 

For  nomenclature,  see  p.  50. 
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Parameters  of  Most  Unstable  Mode 


For  nomenclature,  see  p.  50. 
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In  the  preceding  tables,  subscripts  e,  i,  d,  b,  p  stand  for  “electron,”  “ion,” 
“drift,”  “beam,”  and  “plasma,”  respectively.  Thermal  velocities  are  denoted 
by  a  bar.  In  addition,  the  following  are  used: 


m 

electron  mass 

M 

ion  mass 

V 

velocity 

T 

temperature 

44  g  5  Tl  i 

number  density 

Tl 

harmonic  number 

Cs  =  (Te/ M)1/2 

ion  sound  speed 

^  e  j  ^  i 

plasma  frequency 

A D 

Debye  length 

re,ri 

gyroradius 

P 

plasma/magnetic  energy 

density  ratio 

VA 

Alfven  speed 

gyrofrequency 

Qh 

hybrid  gyrofrequency, 

Qh2  = 

U 

relative  drift  velocity  of 

two  ion  species 


APPROXIMATE  MAGNITUDES 
IN  SOME  TYPICAL  PLASMAS 


Plasma  Type 

_  G 

n  cm 

T  eV 

LOpe  sec  1 

\e>  cm 

n\r>3 

uei  sec  ' 

Interstellar  gas 

i 

1 

6  x  104 

7  X  102 

4  x  10s 

7  x  10-5 

Gaseous  nebula 

103 

1 

2  x  106 

20 

8  x  106 

6  x  10-2 

Solar  Corona 

109 

102 

2  x  109 

2  x  10_1 

8  x  106 

60 

Diffuse  hot  plasma 

1012 

102 

6  x  1010 

7  X  10-3 

4  x  105 

40 

Solar  atmosphere, 

1014 

1 

6  x  1011 

7  X  10-5 

40 

2  x  109 

gas  discharge 

Warm  plasma 

1014 

10 

6  x  1011 

2  x  10-4 

8  x  102 

107 

Hot  plasma 

1014 

102 

6  x  1011 

7  x  10-4 

4  x  104 

4  x  106 

Thermonuclear 

1015 

104 

2  x  1012 

2  x  10-3 

8  x  106 

5  x  104 

plasma 

Theta  pinch 

1016 

102 

6  x  1012 

7  x  10-5 

4  x  103 

3  x  10s 

Dense  hot  plasma 

o 

00 

102 

6  x  1013 

7  x  10-6 

4  x  102 

2  x  1010 

Laser  Plasma 

102° 

102 

6  x  1014 

7  x  10-7 

40 

2  x  1012 

The  diagram  (facing)  gives  comparable  information  in  graphical  form.22 
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LASERS 


System  Parameters 

.  .  o  1 

Efficiencies  and  power  levels  are  approximate. 


Type 

Wavelength 

Efficiency 

Power  levels  available  (W) 

Om) 

Pulsed 

CW 

C02 

10.6 

0.01-0.02 

(pulsed) 

>  2  x  1013 

>  105 

CO 

5 

0.4 

>  109 

>  100 

Holmium 

2.06 

0.03t— 0.1$ 

>  107 

80 

Iodine 

1.315 

0.003 

3  X  1012 

- 

Nd-glass 

1.06 

- 

1.25  X  1015 

- 

Nd:YAG 

1.064 

- 

109 

>  104 

Nd:YLF 

1.045, 

1.54,1.313 

— 

4  x  108 

80 

Nd:YV04 

1.064 

- 

- 

>  20 

Er:YAG 

2.94 

- 

1.5  x  105 

- 

*Color  center 

1-4 

O 

1 

CO 

5  x  108 

1 

*Ti:Sapphire 

0.7-1. 5 

0.4  X  r/p 

1014 

150 

Ruby 

0.6943 

<  10“3 

1010 

1 

He-Ne 

0.6328 

O 

1 

- 

1-50X10-3 

*Argon  ion 

0.45-0.60 

O 

1 

CO 

5  x  104 

150 

*OPO 

0.3-10 

>  0.1  X  rip 

1010 

5 

n2 

0.3371 

0.001-0.05 

106 

- 

*Dye 

0. 3-1.1 

O 

1 

CO 

5  x  107 

>  100 

Kr-F 

0.26 

0.08 

O 

I-1 

to 

500 

Xenon 

0.175 

0.02 

>  10s 

- 

Ytterbium  fiber 

1.05-1.1 

0.55 

5  x  107 

104 

Erbium  fiber 

1.534 

- 

7  x  106 

100 

Semiconductor 

0.375-1.9 

>  0.5 

3  x  109 

>  103 

*Tunable  sources  jlamp-driven  |diode-driven 


Nd  stands  for  Neodymium;  Er  stands  for  Erbium;  Ti  stands  for  Titanium; 
YAG  stands  for  Yttrium— Aluminum  Garnet;  YLF  stands  for  Yttrium  Lithium 
Fluoride;  YV05  stands  for  Yttrium  Vanadate;  OPO  for  Optical  Parametric 
Oscillator;  r)p  is  pump  laser  efficiency. 
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Formulas 


An  e-m  wave  with  k  ||  B  has  an  index  of  refraction  given  by 
n±  =  [1  -  up2e/uj(uj  =p  t oce)]1/2, 

where  ±  refers  to  the  helicity.  The  rate  of  change  of  polarization  angle  0  as  a 
function  of  displacement  s  (Faraday  rotation)  is  given  by 

dO / ds  =  (fc/2)(n_  —  n_|_)  =  2.36  X  10 4 N B f  2  cm  1 , 

where  N  is  the  electron  number  density,  B  is  the  field  strength,  and  /  is  the 
wave  frequency,  all  in  cgs. 

The  quiver  velocity  of  an  electron  in  an  e-m  field  of  angular  frequency  c o 
is 

vo  =  eEmax/ mcj  =  25.6 1  '  Ao  cm  sec 

in  terms  of  the  laser  flux  I  =  cEm 2x/87r,  with  I  in  watt/cm2,  laser  wavelength 
Ao  in  p,m.  The  ratio  of  quiver  energy  to  thermal  energy  is 

Vbqu/Wth  =  mev02 /2kT  =  1.81  X  10_13Ao2//T, 

where  T  is  given  in  eV.  For  example,  if  I  =  1015  Wcm-2,  Ao  =  1  pm,  T  = 
2  keV,  then  VTqu/VTth  ~  0.1. 

Pondermotive  force: 


T  =  NV(E2)/8ttNc, 

where 

Nc  =  1.1  x  1021Ao  2cm  3. 

For  uniform  illumination  of  a  lens  with  /-number  F,  the  diameter  d  at 
focus  (85%  of  the  energy)  and  the  depth  of  focus  l  (distance  to  first  zero  in 
intensity)  are  given  by 

d  «  2A4F\6/6dl  and  l  ^  ±2F2 X6/6DL. 

Here  6  is  the  beam  divergence  containing  85%  of  energy  and  Qdl,  is  the 
diffraction-limited  divergence: 


Odl  =  2.44A/6, 

where  b  is  the  aperture.  These  formulas  are  modified  for  nonuniform  (such  as 
Gaussian)  illumination  of  the  lens  or  for  pathological  laser  profiles. 
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ATOMIC  PHYSICS  AND  RADIATION 


Energies  and  temperatures  are  in  eV;  all  other  units  are  cgs  except  where 
noted.  Z  is  the  charge  state  (Z  =  0  refers  to  a  neutral  atom);  the  subscript  e 
labels  electrons.  N  refers  to  number  density,  n  to  principal  quantum  number. 
Asterisk  superscripts  on  level  population  densities  denote  local  thermodynamic 
equilibrium  (LTE)  values.  Thus  Nn*  is  the  LTE  number  density  of  atoms  (or 
ions)  in  level  n. 

Characteristic  atomic  collision  cross  section: 

(1)  7rao2  =  8.80  X  10  17  cm2. 

Binding  energy  of  outer  electron  in  level  labelled  by  quantum  numbers  n,  l: 

7 jjH 

(2)  E^(n,l)  =  -- - -f-, 

(n  -  A i)2 


where  E ^  =  13.6  eV  is  the  hydrogen  ionization  energy  and  A;  =  0.75 1  5, 
l  5,  is  the  quantum  defect. 

Excitation  and  Decay 

Cross  section  (Bethe  approximation)  for  electron  excitation  by  dipole 
allowed  transition  m  — >  n  (Refs.  32,  33): 


(3) 


where  fmn  is  the  oscillator  strength,  g(n,m )  is  the  Gaunt  factor,  e  is  the 
incident  electron  energy,  and  A Enrn  =  En  —  Em. 

Electron  excitation  rate  averaged  over  Maxwellian  velocity  distribution,  Xmn 
=  Ne((Tmnv)  (Refs.  34,  35): 


(4) 


X 


mn 


1.6  x  10 


—  5  fmn  (g(n,  m))Ne 


AEnmTe 


1/2 


exp 


-l 

sec 


where  (g(n,  m))  denotes  the  thermal  averaged  Gaunt  factor  (generally  ~  1  for 
atoms,  0.2  for  ions). 
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Rate  for  electron  collisional  deexcitation: 


(5)  Ynm  =  (Nm*/ Nn*)Xrnn. 

Here  Nm*  /  Nn*  =  (gm/gn)  exp(ASri,m/Te)  is  the  Boltzmann  relation  for  level 
population  densities,  where  gn  is  the  statistical  weight  of  level  n. 

Rate  for  spontaneous  decay  n  — »  m  (Einstein  A  coefficient)34 

(6)  Anrn  =  4.3  X  10'  (gm / gn)  frnn(AEnrn)2  sec-1. 


Intensity  emitted  per  unit  volume  from  the  transition  n  — >  m  in  an  optically 
thin  plasma: 


(7) 


I  nm  —  1-6  X  10  19  AnrnNnAEnm  watt/cm3. 


Condition  for  steady  state  in  a  corona  model: 

(8)  Nq  N  e  (aonv)  =  NnAn  o, 


where  the  ground  state  is  labelled  by  a  zero  subscript. 

Hence  for  a  transition  n  — >  m  in  ions,  where  ( g(n ,  0))  ~  0.2, 


(9) 


In.m.  =  5.1  X  10 


—  25  f nmdrrt  NeN0  f  A  Er 


9qT( 


1/2 


V  A  E 


exp 


nO 


AEno  \  watt 
Te  )  cm3 


Ionization  and  Recombination 

In  a  general  time-dependent  situation  the  number  density  of  the  charge 
state  Z  satisfies 


(10) 


dN(Z) 

dt 


=  NP 


+S(Z 


-  S{Z)N{Z)  -  a(Z)N (Z) 

1  )N(Z  -  1)  +  a(Z  +  1  )N(Z  +  1) 


Here  S(oZ )  is  the  ionization  rate.  The  recombination  rate  ct(Z)  has  the  form 
ol{Z)  =  ar(Z )  +  N ea.2,(Z) ,  where  ar  and  a 3  are  the  radiative  and  three-body 
recombination  rates,  respectively. 
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Classical  ionization  cross-section36  for  any  atomic  shell  j 


(11) 


cr.i  =  6  X  10  14bjgj(x)/Uj2  cm2. 


Here  bj  is  the  number  of  shell  electrons;  U j  is  the  binding  energy  of  the  ejected 
electron;  x  =  e / Uj ,  where  e  is  the  incident  electron  energy;  and  g  is  a  universal 
function  with  a  minimum  value  gm in  ~  0.2  at  x  ~  4. 

Ionization  from  ion  ground  state,  averaged  over  Maxwellian  electron  distribu¬ 
tion,  for  0.02  <  Te/E ^  <  100  (Ref.  35): 


(12)  S(Z)  =  10 


-5 


(Te/El) 


,1/2 


(^)3/2(6.0  +  Te/^) 


exp  - 


cm 


’/sec, 


where  E^  is  the  ionization  energy. 

Electron-ion  radiative  recombination  rate  (e  +  N (Z)  — »  N (Z  —  1)  +  hu) 
for  Te/Z2  <  400  eV  (Ref.  37): 


(13) 


ar(Z) 


5.2  x  10  14Z 


0.43+  bn^/Te) 


+0.469(£^,/Te) 


sec. 


For  1  eV  <  Te/.Z2  <  15  eV,  this  becomes  approximately35 

(14)  ar(Z )  =  2.7  X  10-13Z2Te-1/2  cm3/sec. 


Collisional  (three-body)  recombination  rate  for  singly  ionized  plasma:38 

(15)  03  =  8.75  X  10  2<Te  4-5  cm6/sec. 


Photoionization  cross  section  for  ions  in  level  n,  l  (short-wavelength  limit): 

(16)  crph(n,Z)  =  1.64  X  10  16Z5/n3if,l2!cm2, 

where  K  is  the  wavenumber  in  Rydbergs  (1  Rydberg  =  1.0974  X  10°  cm-1). 
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Ionization  Equilibrium  Models 

Saha  equilibrium:39 


(17) 


NeN^{Z)  „  afT.3/2 


Nn*(Z  -  1) 


^  n  w  i,2l9l  1e  '  (  E^(n,l) 

=  6.0  X  10  - exp - - - 

9n  V  Je 


cm 


—  3 


where  g ^  is  the  statistical  weight  for  level  n  of  charge  state  Z  and  E^(n,Z) 
is  the  ionization  energy  of  the  neutral  atom  initially  in  level  (n,Z),  given  by 
Eq.  (2). 


In  a  steady  state  at  high  electron  density, 


(18) 


NeN*(Z )  _  S{Z  -  1) 
N*(Z  -  1)  ~~ 


a  function  only  of  T . 

Conditions  for  LTE:39 

(a)  Collisional  and  radiative  excitation  rates  for  a  level  n  must  satisfy 
(19)  Ynm  >  10 Anm. 


(b)  Electron  density  must  satisfy 

(20)  Ne  >  7  X  1018 Z7 n~ 17/2  (T /E^)1/2cm-3. 


Steady  state  condition  in  corona  model: 


(21) 


N  (Z  —  1)  Oir 

N(Z)  ~  S{Z  -  1)  ' 


Corona  model  is  applicable  if40 

(22)  1012t/_1  <  Ne  <  1016Te7/2  cm-3, 

where  tj  is  the  ionization  time. 


56 


Radiation 


N.  B.  Energies  and  temperatures  are  in  eV;  all  other  quantities  are  in 
cgs  units  except  where  noted.  Z  is  the  charge  state  ( Z  =  0  refers  to  a  neutral 
atom);  the  subscript  e  labels  electrons.  N  is  number  density. 

Average  radiative  decay  rate  of  a  state  with  principal  quantum  number  n  is 

(23)  An  =  ^  ^  Anrn  =  1.6  x  10loZ4n  9^2  sec. 

m  <  n 


Natural  linewidth  (A E  in  eV): 

(24)  A E  At  =  h  =  4.14  x  10-15  eV  sec, 

where  At  is  the  lifetime  of  the  line. 

Doppler  width: 

(25)  AA/A  =  7.7  x  10_5(T/p,)1/2, 

where  /j,  is  the  mass  of  the  emitting  atom  or  ion  scaled  by  the  proton  mass. 

on 

Optical  depth  for  a  Doppler-broadened  line: 

(26)  r  =  3.52xlO~13fnmX(Mc2 /kT)1/2NL  =  5.4  x  10~9  fmn  X(/x/T)1/2  NL, 

where  fnm  is  the  absorption  oscillator  strength,  A  is  the  wavelength,  and  L  is 
the  physical  depth  of  the  plasma;  M,  N ,  and  T  are  the  mass,  number  density, 
and  temperature  of  the  absorber;  /u,  is  M  divided  by  the  proton  mass.  Optically 
thin  means  r  <  1. 

Resonance  absorption  cross  section  at  center  of  line: 

(27)  aA=Ac  =  5.6  X  10_13A2/AAcm2. 

Wien  displacement  law  (wavelength  of  maximum  black-body  emission): 

(28)  Amax  =  2.50  X  10-5T-1  cm. 

Radiation  from  the  surface  of  a  black  body  at  temperature  T: 

(29)  W  =  1.03  X  105T4  watt/cm2. 
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Bremsstrahlung  from  hydrogen-like  plasma:26 

(30)  PBr  =  1-69  X  10~32 NeTe1/2  ^  [ Z2N(Z )]  watt/ 


cm 


where  the  sum  is  over  all  ionization  states  Z . 

Bremsstrahlung  optical  depth:41 

(31)  r  =  5.0  X  10~38 NeNiZ2 gLT~7/2 , 

where  ~g  ~  1.2  is  an  average  Gaunt  factor  and  L  is  the  physical  path  length. 

Inverse  bremsstrahlung  absorption  coefficient42  for  radiation  of  angular  fre¬ 
quency  <j u: 

(32)  k  =  3.1  X  10  7  Zne2  In  AT  3^2c o  2(1  —  cj2/cj2)  1^2  cm  1; 

here  A  is  the  electron  thermal  velocity  divided  by  V,  where  V  is  the  larger  of 
lj  and  cjp  multiplied  by  the  larger  of  Ze2 /kT  and  Ti/(mkT)1^2. 

Recombination  (free-bound)  radiation: 


(33) 


Pr  =  1.69  x  10 


-32 


N, 


:Te1/2^^  Z2N(Z)  ^ 


E 


z-i  \  1 


watt /cm" 


Cyclotron  radiation26  in  magnetic  field  B: 

(34)  Pc  =  6.21  x  10  23 B2 NeTe  watt/cm3. 

For  NekTe  =  NikTi  =  B2/16n  ((3=1,  isothermal  plasma),26 

(35)  Pc  =  5.00  X  10  38 N2 T2  watt /cm3. 


Cyclotron  radiation  energy  loss  e-folding  time  for  a  single  electron: 

9.0  X  108B-2 


41 


(36) 


sec, 


2.5  +  7 

2 

where  7  is  the  kinetic  plus  rest  energy  divided  by  the  rest  energy  me  . 
Number  of  cyclotron  harmonics41  trapped  in  a  medium  of  finite  depth  L : 

(37)  mtr  =  (57/3BL)1/6, 


where  (3  =  Sir NkT/B2 . 

Line  radiation  is  given  by  summing  Eq.  (9)  over  all  species  in  the  plasma. 
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ATOMIC  SPECTROSCOPY 


Spectroscopic  notation  combines  observational  and  theoretical  elements. 
Observationally,  spectral  lines  are  grouped  in  series  with  line  spacings  which 
decrease  toward  the  series  limit.  Every  line  can  be  related  theoretically  to  a 
transition  between  two  atomic  states,  each  identified  by  its  quantum  numbers. 

Ionization  levels  are  indicated  by  roman  numerals.  Thus  C  I  is  unionized 
carbon,  C  II  is  singly  ionized,  etc.  The  state  of  a  one-electron  atom  (hydrogen) 
or  ion  (Hell,  Li  III,  etc.)  is  specified  by  identifying  the  principal  quantum 
number  n  =  1,2,...,  the  orbital  angular  momentum  1  =  0,  1,  .  .  .  ,  n  —  1,  and 
the  spin  angular  momentum  s  =  =t^.  The  total  angular  momentum  j  is  the 
magnitude  of  the  vector  sum  of  1  and  s,  j  =  l  =h  \  ( j  >  ^ ).  The  letters  s, 
p,  d,  f,  g,  h,  i,  k,  1,  .  .  .  ,  respectively,  are  associated  with  angular  momenta 
l  =  0,  1,  2,  3,  4,  5,  6,  7,  8,  ...  .  The  atomic  states  of  hydrogen  and  hydrogenic 
ions  are  degenerate:  neglecting  fine  structure,  their  energies  depend  only  on  n 
according  to 

RoohcZ2  n~2  Ry  Z2 

1  +  m/M  n2  ’ 

where  h  is  Planck’s  constant,  c  is  the  velocity  of  light,  m  is  the  electron  mass, 
M  and  Z  are  the  mass  and  charge  state  of  the  nucleus,  and 

Roo  =  109,  737  cm  1 

is  the  Rydberg  constant.  If  En  is  divided  by  he,  the  result  is  in  wavenumber 
units.  The  energy  associated  with  a  transition  m  — >  n  is  given  by 

A Emn  =  Ry(l /m2  -  1/n2), 


with  m  <  n  ( m  >  n )  for  absorption  (emission)  lines. 

For  hydrogen  and  hydrogenic  ions  the  series  of  lines  belonging  to  the 
transitions  m  — »  n  have  conventional  names: 


Transition 

1  — >  n  2  — >  n  3  — »  n  4  — >  n  5  — 6  — >  n 

Name 

Lyman  Balmer  Paschen  Brackett  Pfund  Humphreys 

Successive  lines  in  any  series  are  denoted  a,  (3,  7,  etc.  Thus  the  transition  1  — > 
3  gives  rise  to  the  Lyman -(3  line.  Relativistic  effects,  quantum  electrodynamic 
effects  (e.g.,  the  Lamb  shift),  and  interactions  between  the  nuclear  magnetic 


59 


moment  and  the  magnetic  field  due  to  the  electron  produce  small  shifts  and 

_  Q  _  1 

^  10  cm  ;  these  last  are  called  “hyperfine  structure.” 

In  many-electron  atoms  the  electrons  are  grouped  in  closed  and  open 
shells,  with  spectroscopic  properties  determined  mainly  by  the  outer  shell. 
Shell  energies  depend  primarily  on  n;  the  shells  corresponding  to  n  =  1,  2, 
3,  .  .  .  are  called  K,  L,  M ,  etc.  A  shell  is  made  up  of  subshells  of  different 
angular  momenta,  each  labeled  according  to  the  values  of  n,  l ,  and  the  number 
of  electrons  it  contains  out  of  the  maximum  possible  number,  2(2 1  +  1).  For 
example,  2p5  indicates  that  there  are  5  electrons  in  the  subshell  corresponding 
to  l  =  1  (denoted  by  p)  and  n  =  2. 

In  the  lighter  elements  the  electrons  fill  up  subshells  within  each  shell 
in  the  order  s,  p,  d,  etc.,  and  no  shell  acquires  electrons  until  the  lower  shells 
are  full.  In  the  heavier  elements  this  rule  does  not  always  hold.  But  if  a 
particular  subshell  is  filled  in  a  noble  gas,  then  the  same  subshell  is  filled  in 
the  atoms  of  all  elements  that  come  later  in  the  periodic  table.  The  ground 
state  configurations  of  the  noble  gases  are  as  follows: 

He  Is2 

Ne  ls22s22p6 

Ar  ls22s22p63s23p6 

Kr  ls22s22p63s23p63d104s24p6 

Xe  ls22s22p63s23p63d104s24p64d105s25p6 

Rn  ls22s22p63s23p63d104s24p64d104f  145s25p65d106s26p6 

Alkali  metals  (Li,  Na,  K,  etc.)  resemble  hydrogen;  their  transitions  are  de¬ 
scribed  by  giving  n  and  l  in  the  initial  and  final  states  for  the  single  outer 
(valence)  electron. 

For  general  transitions  in  most  atoms  the  atomic  states  are  specified  in 
terms  of  the  parity  (  —  l)^li  and  the  magnitudes  of  the  orbital  angular  momen¬ 
tum  L  =  El,-  ,  the  spin  S  =  Es^,  and  the  total  angular  momentum  J  =  L  +  S, 
where  all  sums  are  carried  out  over  the  unfilled  subshells  (the  filled  ones  sum 
to  zero).  If  a  magnetic  field  is  present  the  projections  Ml,  Ms,  and  M  of 
L,  S,  and  J  along  the  field  are  also  needed.  The  quantum  numbers  satisfy 
|Ml|  <  L  <  ul,  |  Mg  |  <  S  <  i//2,  and  \M\  <  J  <  L  +  S,  where  v  is  the 
number  of  electrons  in  the  unfilled  subshell.  Upper-case  letters  S,  P,  D,  etc., 
stand  for  L  =  0,  1,  2,  etc.,  in  analogy  with  the  notation  for  a  single  electron. 
For  example,  the  ground  state  of  Cl  is  described  by  3pJ  2P3/2-  The  first  part 
indicates  that  there  are  5  electrons  in  the  subshell  corresponding  to  n  =  3  and 
1  =  1.  (The  closed  inner  subshells  Is2 2s2 2p6 3s2,  identical  with  the  configura¬ 
tion  of  Mg,  are  usually  omitted.)  The  symbol  ‘P’  indicates  that  the  angular 
momenta  of  the  outer  electrons  combine  to  give  L  =  1.  The  prefix  ‘2’  repre¬ 
sents  the  value  of  the  multiplicity  2 S  +  1  (the  number  of  states  with  nearly  the 
same  energy),  which  is  equivalent  to  specifying  S  =  The  subscript  3/2  is 
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the  value  of  J .  The  superscript  ‘o’  indicates  that  the  state  has  odd  parity;  it 
would  be  omitted  if  the  state  were  even. 

The  notation  for  excited  states  is  similar.  For  example,  helium  has  a  state 
ls2s  3Si  which  lies  19.72  eV  (159,856  cm-1)  above  the  ground  state  Is2  1So- 
But  the  two  “terms”  do  not  “combine”  (transitions  between  them  do  not  occur) 
because  this  would  violate,  e.g.,  the  quantum-mechanical  selection  rule  that 
the  parity  must  change  from  odd  to  even  or  from  even  to  odd.  For  electric 
dipole  transitions  (the  only  ones  possible  in  the  long- wavelength  limit),  other 
selection  rules  are  that  the  value  of  l  of  only  one  electron  can  change,  and  only 
by  A l  =  dhl;  A S  =  0;  A L  =  ±1  or  0;  and  A  J  =  ±1  or  0  (but  L  =  0  does  not 
combine  with  L  =  0  and  J  =  0  does  not  combine  with  J  =  0).  Transitions 
are  possible  between  the  helium  ground  state  (which  has  S  =  0,  L  =  0,  J  =  0, 
and  even  parity)  and,  e.g.,  the  state  ls2p  1P))  (with  S  =  0,  L  =  1,  J  =  1, 
odd  parity,  excitation  energy  21.22  eV).  These  rules  hold  accurately  only  for 
light  atoms  in  the  absence  of  strong  electric  or  magnetic  fields.  Transitions 
that  obey  the  selection  rules  are  called  “allowed”;  those  that  do  not  are  called 
“forbidden.” 

The  amount  of  information  needed  to  adequately  characterize  a  state  in- 
creases  with  the  number  of  electrons;  this  is  reflected  in  the  notation.  Thus 
O II  has  an  allowed  transition  between  the  states  2p  3p 
2 -^7/2  an<^  2p2(1D)3d/  2F7/2  (and  between  the  states  obtained  by  changing 
J  from  7/2  to  5/2  in  either  or  both  terms).  Here  both  states  have  two  elec¬ 
trons  with  n  =  2  and  l  =  1;  the  closed  subshells  ls22s2  are  not  shown.  The 
outer  (n  =  3)  electron  has  l  =  1  in  the  first  state  and  l  =  2  in  the  second. 
The  prime  indicates  that  if  the  outermost  electron  were  removed  by  ionization, 
the  resulting  ion  would  not  be  in  its  lowest  energy  state.  The  expression  (1D) 
give  the  multiplicity  and  total  angular  momentum  of  the  “parent”  term,  i.e., 
the  subshell  immediately  below  the  valence  subshell;  this  is  understood  to  be 
the  same  in  both  states.  (Grandparents,  etc.,  sometimes  have  to  be  specified 
in  heavier  atoms  and  ions.)  Another  example43  is  the  allowed  transition  from 
2p2(3P)3p  2Pi/2  (or  2P|j/2 )  to  2p2(1D)3d/  in  which  there  is  a  “spin 

flip”  (from  antiparallel  to  parallel)  in  the  n  =  2,  l  =  1  subshell,  as  well  as 
changes  from  one  state  to  the  other  in  the  value  of  l  for  the  valence  electron 
and  in  L. 

The  description  of  fine  structure,  Stark  and  Zeeman  effects,  spectra  of 
highly  ionized  or  heavy  atoms,  etc.,  is  more  complicated.  The  most  important 
difference  between  optical  and  X-ray  spectra  is  that  the  latter  involve  energy 
changes  of  the  inner  electrons  rather  than  the  outer  ones;  often  several  electrons 
participate. 
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COMPLEX  (DUSTY)  PLASMAS 


Complex  (dusty)  plasmas  (CDPs)  may  be  regarded  as  a  new  and  unusual 
state  of  matter.  CDPs  contain  charged  microparticles  (dust  grains)  in  addition 
to  electrons,  ions,  and  neutral  gas.  Electrostatic  coupling  between  the  grains 
can  vary  over  a  wide  range,  so  that  the  states  of  CDPs  can  change  from  weakly 
coupled  (gaseous)  to  crystalline.  CDPs  can  be  investigated  at  the  kinetic  level 
(individual  particles  are  easily  visualized  and  relevant  time  scales  are  accessi¬ 
ble).  CDPs  are  of  interest  as  a  non-Hamiltonian  system  of  interacting  particles 
and  as  a  means  to  study  generic  fundamental  physics  of  self-organization,  pat¬ 
tern  formation,  phase  transitions,  and  scaling.  Their  discovery  has  therefore 
opened  new  ways  of  precision  investigations  in  many-particle  physics. 

Typical  experimental  dust  properties 

grain  size  (radius)  a  ~  0.3  —  30  pm,  mass  md  ~  3  X  10“  '  —  3  X  10“  13  g,  number 
density  (in  terms  of  the  interparticle  distance)  rid  103  -  107  cm-3, 

temperature  Td  ~  3  X  10-2  —  102  eV. 

Typical  discharge  (bulk)  plasmas 

gas  pressure  p  ~  10"  2  —  1  Torr,  Tj  ~  Tn  ~  3  X  10“  2  eV,  VTi  —  7  X  104  cm/s 
(Ar),  Te  ~  0.3  —  3  eV,  rii  ~  ne  ~  108  —  1010  cm-3,  screening  length  Ad  ~ 
A r>i  ~  20  —  200  fi m,  u>pi  ~  2  x  106  —  2  x  107  s_1  (Ar).  B  fields  up  to  B  ~  3  T. 

Dimensionless 


Havnes  parameter 

P  = 

\Z\rid/ne 

normalized  charge 

z  = 

\Z\e2  /  kTea 

dust-dust  scattering  parameter 

Pd  = 

=  Z2e2/kTd\D 

dust-plasma  scattering  parameter 

Pe,i 

=  \Z\e2 /kTejiX]j 

coupling  parameter 

r  = 

(Z2e2 /kTdA)  exp(-A/Ac) 

lattice  parameter 

K  = 

A/ 

particle  parameter 

a  = 

a/ A 

lattice  magnetization  parameter 

p  = 

A /rd 

Typical  experimental  values:  P  10-4  —  102 ,z  ~  2-4  (Z  ~  103  -  105  electron 
charges),  P  <  103,  n  ~  0.3  —  10,  a  ~  10-4  —  3  X  10-2,  p  <  1 


Frequencies 

dust  plasma  frequency 


charge  fluctuation  frequency 


Upd  =  (47rZ2e2nd/md)1/2 
—  {\Z\TPpmi/md)1/2id 

^ch  —  7/2tt  (a/^-D  )^pi 
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dust-gas  friction  rate 
dust  gyrofrequency 


und  ~  10 a2p/mdvTn 
ujcd  =  ZeB/mdc 


Velocities 

dust  thermal  velocity 

dust  acoustic  wave  velocity 

dust  Alfven  wave  velocity 
dust-acoustic  Mach  number 
dust  magnetic  Mach  number 

dust  lattice  (acoustic)  wave  velocity 


VTd  =  (kTd/mdf^  = 

C  da  =  Upd^D 

-  (\Z\-^^mi/md)1/2vTi 
VAd  =  B  /  (Airndmd)1/2 
V/Cda 
V/vAd 

Cql  =  Upd*DFi)t(K) 


The  range  of  the  dust-lattice  wavenumbers  is  KA  <  n  The  functions  F^^k) 
for  longitudinal  and  transverse  waves  can  be  approximated44,45  with  accuracy 
<  1%  in  the  range  k  <  5: 

Fi  ~  2.70k1/2(1  -  0.096k  -  0.004k2),  Ft  ~  0.51k(1  -  0.039k2), 


Lengths 

frictional  dissipation  length  =  VTd/vnd 

dust  Coulomb  radius  Rce,i  =  \Z\ e2 /kTe>i 

dust  gyroradius  rd  =  vTd/ucd 

Grain  Charging 

The  charge  evolution  equation  is  d\Z\/dt  =  Ii  —  Ie.  From  orbital  motion 
limited  (OML)  theory46  in  the  collisionless  limit  ien(in)  ^  cl: 


Ie 


Ii 


\/  87Ti 


7TCL  rtiVT- 


Grains  are  charged  negatively.  The  grain  charge  can  vary  in  response  to  spatial 
and  temporal  variations  of  the  plasma.  Charge  fluctuations  are  always  present, 
with  frequency  cnch-  Other  charging  mechanisms  are  photoemission,  secondary 
emission,  thermionic  emission,  field  emission,  etc.  Charged  dust  grains  change 
the  plasma  composition,  keeping  quasineutrality.  A  measure  of  this  is  the 
Havnes  parameter  P  =  \Z\nd/ne.  The  balance  of  Ie  and  A  yields 

exp(“2)=  teS)  7  (1  +  ^")[1  +  p(z)] 
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When  the  relative  charge  density  of  dust  is  large,  P  >  1,  the  grain  charge  Z 
monotonically  decreases. 

Forces  and  momentum  transfer 

In  addition  to  the  usual  electromagnetic  forces,  grains  in  complex  plasmas  are 
also  subject  to:  gravity  force  Fg  =  m^g;  thermophoretic  force 

4\/27r  2  z  x 

Fth  = - -  (a  /vTn)KnVTn 

15 

(where  Kn  is  the  coefficient  of  gas  thermal  conductivity);  forces  associated 
with  the  momentum  transfer  from  other  species,  Fa  =  -m^^adVad,  i.e., 
neutral,  ion,  and  electron  drag.  For  collisions  between  charged  particles,  two 
limiting  cases  are  distinguished  by  the  magnitude  of  the  scattering  parameter 
(3a.  When  (3a  1  the  result  is  independent  of  the  sign  of  the  potential.  When 

/ 3a  1,  the  results  for  repulsive  and  attractive  interaction  potentials  are 

different.  For  typical  complex  plasmas  the  hierarchy  of  scattering  parameters 
is  (3e( ~  0.01  —  0.3)  <C  1  —  30)  <C  Pd(^  103  —  3  X  104).  The  generic 

expressions  for  different  types  of  collisions  are47 

vad  =  (4\/27t/3  ){m0,/md)a2n0iVTa^ad 

Electron-dust  collisions 


1  2 

4>ed  ~  -Z  Aed  Pe  <C  1 

Zj 

Ion-dust  collisions 

$  =  f  iz2(Te/Ti)2Aid  Pi  <5 

zd  {2(An /a)2  (In2  Pi +2  In  Pi +  2),  Pi  >  13 

Dust-dust  collisons 

Q  —  /  Zd^dd  Pd  ^  1 

\  (AD/a)2  [In  4pd  -  In  In  4/3d] ,  Pd  >  1 


where  zd  =  Z2e2/akTd- 


For  t>dd  ~  vnd  the  complex  plasma  is  in  a  two-phase  state,  and  for  und  !$>  Vdd 
we  have  merely  tracer  particles  (dust-neutral  gas  interaction  dominates).  The 
momentum  transfer  cross  section  is  proportional  to  the  Coulomb  logarithm 
A ad  when  the  Coulomb  scattering  theory  is  applicable.  It  is  determined  by 
integration  over  the  impact  parameters,  from  pm in  to  pm ax-  Pmin  is  due  to  finite 
grain  size  and  is  given  by  OML  theory.  pm ax  =  for  repulsive  interaction 
(applicable  for  pa  -C  1),  and  pm ax  =  An  (1  +  2poi  )1//2  for  attractive  interaction 
(applicable  up  to  Pa  <  5). 
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For  repulsive  interaction  (electron-dust  and  dust-dust) 


A-atd 


ln[l  +4(A e>  / a a)2 x2]dx 


2  zf 


1  )dx, 


where  ze  =  z,  ae  =  a,  and  =  2 a. 
For  ion-dust  (attraction) 


A-id 


1  +  2{Ti/Te){\D /a)x 
1  +  2 (Ti /Te)x 


dx. 


For  Udd  v-nd  the  complex  plasma  behaves  like  a  one  phase  system  (dust-dust 
interaction  dominates). 

Phase  Diagram  of  Complex  Plasmas 

The  figure  below  represents  different  “phase  states”  of  CDPs  as  functions  of 
the  electrostatic  coupling  parameter  F  and  re  or  a,  respectively.  The  verti¬ 
cal  dashed  line  at  re  =  1  conditionally  divides  the  system  into  Coulomb  and 
Yukawa  parts.  With  respect  to  the  usual  plasma  phase,  in  the  diagram  be¬ 
low  the  complex  plasmas  are  “located”  mostly  in  the  strong  coupling  regime 
(equivalent  to  the  top  left  corner). 

Regions  I  (V)  represent  Coulomb  (Yukawa)  crystals,  the  crystallization  condi¬ 
tion  is48  r  >  106(1  +  re  +  re2/2)_  1  .  Regions  II  (VI)  are  for  Coulomb  (Yukawa) 
non-ideal  plasmas  —  the  characteristic  range  of  dust-dust  interaction  (in  terms 
of  the  momentum  transfer)  is  larger  than  the  intergrain  distance  (in  terms  of 
the  Wigner-Seitz  radius),  (cr/ 7r)1//2  >  (47r/3)-1//3A,  which  implies  that  the 
interaction  is  essentially  multiparticle. 

-i  .  , 

Regions  III  (VII  and  VIII)  correspond  to 
Coulomb  (Yukawa)  ideal  gases.  The  range 
of  dust-dust  interaction  is  smaller  than  the 
intergrain  distance  and  only  pair  collisions 
are  important.  In  addition,  in  the  region 
VIII  the  pair  Yukawa  interaction  asymp¬ 
totically  reduces  to  the  hard  sphere  limit, 
forming  a  “Yukawa  granular  medium” .  In 
region  IV  the  electrostatic  interaction  is 
unimportant  and  the  system  is  like  a  uaual 
granular  medium. 


a  =/ 
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